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DEVELOR€ENT OF A COMPUTER CODE FOR CALCULATING THE 
STEADY SUPER/HYPFRSONIC IHVISCID FLOW AROUND 
— REAL CONFIGURATIONS- - — — 


VOLUME II - CODE DESCRIPTION 

by 

F. Marconi and L.S. Yaeger 
GRUMMAN AEROSPACE CORPORATION 


SUMMARY 

A set of four computer codes has been developed to compute the invis- 
cid super/hypersonic flow field about complex vehicle geometries. The 
numerical procedures used in these codes are described in detail in Volme 
I of this report. Here the codes developed are described with two views; 
one oriented toward the user and the other toward the programmer. 

The nomenclature used in the codes, the input and output formats, and 
the storage requirements and computer time are discussed in detail. A 
description of routines, over-all logic flow, and overlay structure are 
also presented. 
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ORIENTATION 


Volume I gave the approach to the numerics and this volume gives all 
the other matters closely related to the codes. 

When handling real configurations, one matter of decisive importance 
is how the geometiy is modeled, particularly when one solves partial dif- 
ferential equations (rather than integral equations as in other aero- 
dynamic efforts ) . Therefore the user needs to have an idea of the approach 
to the geometry modeling before learning the operations of the codes. In 
this piece of work, geometry modeling is done with a technique developed 
by A. Vachris and L. Yaeger. For the reader unfamiliar with it. Appendix A 
gives a brief, self-contained description of this technique, called the 
QUICK Geometry System. Appendix A is couched in code oriented terms with- 
out indulging in dissertations of lofting techniques. 

To compute the transonic flow over the nose of blunted vehicles, a 
three dimensional time asymptotic technique was used. The code (BLUNT) 
which was used for these calculations is briefly discussed in Appendix B. 
The computational procedxire used is discussed in reference 1. 

The typical user will be interested in Part 1, the 'user oriented... 
documentation', which will give him the minimm amoiuit of . information 
necessary to operate, as 'black boxes', the codes developed or adapted 
under this contract. The large amount of nomenclatTxre included here is 
to be used primaj’ily as a dictionary; only a few symbols and terms need to 
be learned by the user, namely those that appear in the input/output data 
format. The nomenclature and non-dimensionalizations of Volume I of this 
report are used here. 

To the programmer who wants to look into the 'black boxes'. Part 2 
of this volume is dedicated. 



PART 1 USER-ORIENTED DOCUMENTATION 
DESCRIPTION OF CODES AND THEIR INTERACTION 

A series of five codes has been developed or adapted under this 
contract: 

QUICK written by A. Vachris and L.S. Yaeger, is a geometry system 

designed to allow the user to model a complex vehicle geo- 
metry in a quick, straightforward fashion. The QUICK geo- 
metry system also allows another code, which uses the 
modeled vehicle geometry as input, to interrogate the model 
for cross sectional Information as efficiently as possible. 
QUICK consists of an initial defining and logical checkout 
group of routines , which actually set up the mathematical 
model, and a second group of routines (called SUB-QUICK 
throughout this report) which is used for interrogating the 
mathematical model. SUB-QUICK is used as a palrt of QUICK 
to inspedt the modeled vehicle, and as a part of the supersonic 
flow field code (STEIN), along with an output data set frcM 
QUICK (the QUICK intermediate data deck), to supply all 
geometry information. 

STEIN - written by F. Marconi and L.S. Yaeger, is a supersonic flow 
field code designed as a tool to allow the user to compute 
the super/hypersonic inviscid flow about realistic configura- 
tions. The numerical techniques utilized in STEIN are 
described in detail in Volume I of this report. STEIN reads 
control data, starting plane data, and geometry data, and 
computes the flow from the starting plane to a user prescribed 
axial station. The nose region of the vehicle must be com- 
^ puted with another code which generates starting plane data 

(where the axial Mach number is supersonic, see Volume l). 

In STEIN there is a routine which will compute the starting 
plane data for sharp circular cones at small angles of attack, 
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so the initial data need not he • generated elsewhere for this 
case. For blunt nose vehicles a BLUNT BODY code, developed 
by Professor Gino Moretti (ref. 1) which is compatible with 
~ -QUICK and STEIN is used. STEIN compute^ the flow field, the 

aerodynamic coefficients and the metric coefficient from 
the starting plane to the end station. 

STIIMBL - written by L.S. Yaeger, is a code designed to utilize flow 
field data, output on tape from STEIN, to compute stream- 
lines on the body, create pseudostream surfaces (p-s-s; 
defined by the body surface normals taken at each point 
along a given body streamline), and evaluate flow variables 
and their normal derivatives along the streamlines and in 
' the p-s-s from the starting plane to the end station. 

BOOM - written by L.S. Yaeger, is a code designed to utilize flow 

field data (from the same STEIN output tape used by STRMBL) 
to evaluate flow variables on a data cylinder (whose center- 
line is the z-axis and radius is user-specified) for sonic 
boom work. 

BLUNT - developed by Moretti, uses a time dependent computational 

technique to asymptote to a steady transonic solution. 

Its results are used as an initial condition to compute 
three dimensional supersonic flow over blunt nose vehicles. 
Details of the technique used to compute the blunt nose 
flow fields are presented in reference 1. The geometry 
input for this blunt body code can be either supplied by 
the geometry package ("QUICK") or computed internally for 
simple noses. The output from this code is compatible with 
the three dimensional supersonic flow field code's (STEIN) 
requirements for initial data. The input for BLUNT is 
described in Appendix B. 

The interaction of these codes (i.e., input-output flow) is described in 

figure 1. 
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Figure 1 - INTERACTION OF SYSTEM OF CODES 
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NOMENCIA'rURE 


QUICK TERMINOLOGY 

During the discussion of QUICK, several terms will appear frequently, 
and as such, will he defined here: 

1) Cross section - standard definition; a planar cut through the 
vehicle normal to the FRL at a given x- station. 

2) Cross-sectional model - mathematical abstraction of a cross 
section, using simple curves to represent arcs between specified 
control points. 

3) Control points - break or joining points for defining each arc. 

4) Arc - a portion of one simple mathematical curve between two 
control points in cross section. 

5) Body lines - the defining lines of. the vehicle geometry in plan 
and profile views; x-running control points given as y^ = y^(x) 
and/or z.=z.(x). 

6) Body line model - mathematical abstraction of a body line, using 
simple curves to represent se@nents between specified match points 

7) Match or Key points - break or joining points between body line 
segments; initial and terminal points for defining each segment. 

8) Segment - a portion of one simple mathematical curve between two 
match points of a body line model. 

9) Component - same as an arc; usually considered to be a named 
portion of the vehicle geometry (e.g., a wing-upper- ellipse may 
be component WNGUPELL). 

Body line segments are discussed in terras of an origin point at (x^, v^) 

(v standing for y or z), a termination point (xg, Vg), an initial slope t^ 
and a final slope t^. 
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SYMBOL LIST FOR QUICK 


ANAME 




ARCNM(l) 


ARCNM(2) 


ASHAPE 

ASPEC(l) 


Hollerith input variable; body line (BL) /control point name 
to which BNAME is to be aliased, -when applicable (blank when 
not ) 

Hollerith input variable; cross section (CS) arc or com- 
ponent name 

Hollerith input variable; if type is FILET: the name of the 

most' aft component arc to which the current arc's forward 
end is to be filleted 

If type is other: the name of the most aft component arc 

which, in case of intersection with the current arc, is to 
update the forward end of the current arc and the aft end 
of the intersected arc 

Hollerith input variable; if type is FILET: the name of 

the most forward component arc to which the current arc's 
aft end is to be filleted 

If type is other: the name of the most forward component 

arc which, in case of intersection with the current arc, 
is to update the aft end of the current arc and the forward 
end of the intersected arc 

Hollerith input variable; arc or component shape 
Hollerith input variable: 

= blank yields no effect 

■v 

f • 

= Y when type is FILET, and only y-values are to be specified 
for the next control point in order of input (z is com- 
puted on controlling component) 

= Z when type is FILET, and only z-values are to be specified 
for the next control point in order of input (y is com- 
puted on controlling component) 
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ASPEC(2) 

ATYPE 

AYORZ 

BLCOEF(l 


BLMDEE(I 


= B to indicate that the next control point is the bottom 
centerline of the vehicle for the model cui'rently being 
defined (optional) 

= T to indicate that the next control point is the top 
centerline of the vehicle for the model currently being 
defined (optional) 

Same as ASPEC(l) 


Hollerith input variable; arc or ccmponent type 


Hollerith input variable; the letter Y or Z to indicate 
which definition is to be used when aliasing (blank when 


not) 


1 = 1 

to 7; 

ment and BL 

1 = 1: 

^1 

I = 2: 

^1 

I = 3: 

a2 

1 = 4: 

b" 

1 = 5: 

c 

1 = 6: 

^2 

(H 

II 

^ 2 . 

1 = 1 

to 8; 

1 = 1: 

^1 

I = 2: 

^1 

H 

II 

U) 

X2 

II 

M 


1 = 5: 

^3L 

1 = 6: 

^3L 
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BLMMAX(l) 

BLMMIN(I) 

blmnam(m) 

BLMNYZ(M) 

BNAME 

BTITLE(I, 

BYORZ 

COMPNM(I) 

CPNTNM(I) 

CTITLE(I, 

D(l) 

D(2) 


II 

M 

^3R 

H 

II 

OO 

^3R 


(Xi, v^) and (x^j are the initial and final points, 
respectively, of the given segment, , v ) establishes 

Jij 0-L', 

the slope at the initial side, (x_ , v ) establishes the 

oh 

slope at the terminal side 


I = 1 to KNTBLM; maximum x for each BL model 
I = 1 to KNTBLM; minimum x for each BL model 


Alphanumeric name of each BL model 


Alphanumeric y or z coordinate specification for each BL 
model 

Hollerith input variable; body line/control point name 
which is to be defined 


II ) not used currently 

Hollerith input variable; the letter Y or Z to indicate 
which data coordinate definition is to follow 

I = 1 to KCOMP; component names (alphanumeric) 

I = 1 to KCPNT; control point names (alphanumeric) 

K) I = 1 to 10; alphanumeric CS model title or comments 

Input variable; if type is PIECE or FLINK, this is x^. 
If type is ALINK, PATCH, or FILET, this is a floating 
point mamber equal to KSEG of the segment from which 
x^ and/or v^ are to be determined. 

Input variable; if type is PIECE or FLINK, this is v^. 
If type is ALINK, PATCH, or FILET, this is a floating 

i 

point number equal to KSEG of the segment from which t^ 
is to be determined. 
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\ 

^ D(3) 


D(4). 


D(5) 


B(6) 


HDEL 

HEND 

HGO 

KNOW 

HNOWR 

lAMD 

lANDV 


Input variable; if type is PIECE or ALINK, this is x^. 

If type is FLINK,- PATCH, or FILET, this is a floating point 
number equal to KSEG of the segment from which Xg and/or 
are to be determined. 

Input variable; if type is PIECE or ALINK, this is v^. 

If type is FLINK, PATCH, or FILET, this is a floating point 
number equal to KSEG of the segment. from which is to be 
determined. 

Input variable; if SLPl is blank: 

• -.4 r . 

If type is FILET, this is y^ and t^^ are to be determined 
from the segment specified by D(l) and D(2). If type is 
other, this is x^- . ■ 

If SLPl is other than blank, see definition of SLPl. 

Input vau’iable; if SLP2 is blank: ’ ' 

If type is FILET, this is x^; a>nd t^ are to be determined 
from the segment specified by D(3) and D(4). 

If type is other, this is x^. ■ 

If SLP2 is other than blank, see definition of SLP2. 

Input variable; increment size in degrees to establish inter- 
rogation points between HGO and HEND; not required for 
modes 1 or 3» 

Input variable; final value of theta (in degrees) to be 
interrogated; not required for modes 1 or 3« 

Input variable; initial value of theta (in degrees) to be 
interrogated; not required for modes 1 or 3* 

Current value of 0 in degrees (used in various exercising 
routines; e.g., MODEl, M0DE2, etc). 

Current value of 9 in radians (used in various exercising 
routines, e.g., MODEl, M0DE2, etc). 

lABS(MODE) 

lABS(NDERV) 


9 



1BLC0R(I,J) 

IBLMIX(M) 

IBLMWD(I,W,M) 

IBLMX(I) 

IBLSSH(N,M) 

IBLSX(I) 

IC0MPX(J,K) 

ICRITE 

ICSACC(I,J,K) 


ICSACP(I,J,K) 


/ 

I = 1 to 6 ; body line coordinate index for Yl(l = l), ^ 

Z1(I = 2), Y2(I =3), Z2(I = 4), Y3 and/or Y4(I = 5), Z3 

and/or Z4(l = 6). 

Index to the control point coordinate for which this BL 
model was first defined. 

I = 1 to 4; indicator for the shape (l = l), type (l =2), 
mode of definition (l = 3)» anti freed constraints (l = 4) 
of each segment and BL model 

4* Vi 

I = 1 to NBLCOR; index of the body line model for the I 
coordinate control point 

Shape index for each segment and BL model 

= (1) LINE, (2) CIRC - not used, (3) ELLX, (4) ELLY, 

(5) XPAR, (6) YEAR, (7) RXPA, (8) RYPA, (9) CUE!, (lO) ALL - 
not used, (U) NULL 

I = 1 to KNTBLM: current segment number index for each BL 
model. 

Index of the component definition for each arc and CS model 

Output unit for error and checking messages, primarily for 
use on a time sharing computing system, otherwise, ICRITE 
= IRITE 

I = 1,2; controlling component index for each arc and CS 
model 

1=1: information pertains to forward end of arc 

1=2: information pertains to aft end of arc 

= -1: end of arc is unaffected 

> 0: gives index of another arc which is to intersect the 

th 

J arc for growing pieces, or which is to supply filleting 
th 

information if J arc is a fillet. 

I = 1 to 3: control point index for each arc 6ind CS model 

1=1: initial point of arc 

1=2: final point of arc 

1=3: slope control point for arc 
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\ ICSAFR(J,K) Free constraint index for each arc and CS model (not 
' currently used) 

ICSASH(J,K) Shape index for each arc and CS model 

_ , ^ ^ ^ - LINE - - -* ^ - -- -- - ^ - - 

= 2: CIRC, circle (not currently available) 

=3: ELLI, in-facing ellipse (concave to origin) 

= 4: ELLO, out-facing ellipse (convex to origin) 

ICSASQ(J,K) Sequencing index to establish the order in which cross 
sectional arcs are to be defined. 

ICSATY(J,K) Type index for each arc and CS model 

= 1: PIEC, piece 

= 2: FUN, forward link 

=3: ALIN, aft link 

= 4: PATC, patch 

= 5: FILE, fillet 

ICSMX(KMODEL) Index of current CS model (from 1 to NCSM), describes use 
of library of CS models as applied to this vehicle. 

IFREE Input variable; index of the datum quantity which is to be 

"free," i.e., determined by the code. IFREE ranges from 
1 to 6 corresponding to x^, v^, x^, v^, t^, t^, as ordered. 

A line must have any one of these free; an x- or y- parabola 
must have either 5 or 6 free; other curves should have 
IFREE = 0. 

IN(J ) Indicator^ for each, .arc of , the current CS model 

= -1; arc not included at this station 
= 1: arc included at this station 

IPLOT l/O unit for plot mode output from GEMCHK, MODEl, M0DE2, etc. 

IPUNCH Assumed punch unit (= 7) 


IREAD Input unit 


IRITE Output unit 
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ISPEC(I,J,K) 


ITAPE 

IUOmDR( J ) 

IZBDEX(K) 

IZCDEX 

IZTDEX(K) 

J 

JSEQ 

K 

KARC 

KCOMP 

KCPNT 

KDUM 


I = 1, 2; index to indicate what coordinate is to be speci- 
fied at the initial control point (l = l) and the final 
control point (l = 2) 

=1; y is to be specified (z is to be computed on the 
controlling component) 

=2: z is to be specified (y is to be computed on the 

controlling component)- 

= -1: for nonfillets . 

l/O unit for QUICK intermediate data deck (math model) 

(note: called INREAD in GEOMIN) 

Use order index to establish sequence of CS arcs after inter 
sections and filets are completed. 

Index of the bottom center body line model for each CS model 

Index of the center body line model (mapaxis) 

Index, of the top center body line model for each CS model. 

Index of current cross sectional arc for a given CS model 
(K) from 1 to KWTCSA(K). 

Input variable; definition sequence (order in which the CS 
arcs are to be defined) 

Index of current cross sectional definition (library) model 
(from 1 to NCSM) 

Input variable; number of arcs in current cross sectional 
model. 

** ■ *.) r J i ■ 1 i;,' , 

Number of components used to define all CS models (entire 
vehicle) . 

Number of control points used to define all CS models 
(entire vehicle). 

Input variable; running count of the current cross section 
model . 
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KMODEL 

KNTARC 

KNTBLM 

KNTBLS(M) 

KNTCSA(K) 

KNTCSM 

KSEG 


KZBDEX 

KZCDEX 

KZTDEX 

M 

MODE 


Index of current cross sectional use model (from 1 to 
KWTCSM) 

Number of arcs in the CS model corresponding to the current 

— ' 'Station"" ^ """ " ~ - "" " ' ’ ^ ' 'T 

Number of body line models 

Number of segments for each body line model 

Number of arcs for each cross sectional model 

Number of applications of- cross section models to define 
entire vehicle 

Input variable; the order (in increasing x) in which this 
segment appears in this body line model. A KSEG = -1 
(further arguments not required) terminates the data for 
a given body line. 

Control point index for bottom centerline. 

Control point index for mapaxis. 

Control point index for top centerline, 

Index of current body line definition model (from 1 to 
KNTBLM) 

Input variable; 

= +1, creates body line traces 
= +2, creates cross sectional cuts 

= +3 3 interrogates cross sections in neighborhood of 
control points 

-3 3 allows multiple body line traces to create plan 
and profile views 

= +4, comparison of analytic derivatives with numerically 
formed derivatives 

= +53 check of unit vectors normal to body surface 
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MODEL 

N 

NBLCOR 

NCSM 

NDERV 


NHPTS 

NXPTS 

PNTNAM(l) 


= +6, exercises modes 1, 2, and 3 at the limits of each 
cross sectional model 

-6, exercises modes -2 and -7 at the limits of each cross 
sectional model 

= -7, (plotting mode only) creates cross sectionaJ. cuts, 
but includes all arcs in their entirety (including growing 
pieces still contained within the basic skin) 

Index to the current CS library model definition 

Index of current body line segment for a given BL model (m) 
from 1 to KNTBLS(M). 

Number of control point coordinates to define entire vehicle 
(y and z are distinct, thus NBLCOR = 2*KCPNT). 

Input variable; number of distinct cross section models. 

Input variable,; 

= +N, where N is the order of derivative to be calculated 
(N = 0, 1, or 2) 

= +N, should always be used for checkout interrogations (means 
. each call to a given location is new, thus the radius and 
all temporary variables must be computed) 

= -N, should not be used for checkout interrogations; requires 
previous call to same location (x and 0); radius and certain 
temporary variables are not recomputed. 

Number of Q points ‘(used in' various exercising routines; 
e.g., MODEl, M0DE2, etc). 

Number of x-stations (used in various exercising routines; 
e.g., MODEl, M0DE2, etc). 

Hollerith input variable; control point name for the begin- 
ning of the arc currently being defined. 



PNTNAM(2) 

PNTNAM(3) 

SDEF 

SLPl 


SLP2 


S SHAPE 

STYPE 

THETAl(j) 

THETA2(j) 

TITLE 

UNX 

UNY 

UNZ 


Hollerith input variable; control point name for the termi- 
nation of the arc currently being defined. 

Hollerith input variable; slope control point name for 
the current arc when required, blank if not. 

Hollerith input variable; segment definition mode (currently, 
only two point, two slope/slope control point method is 
available - input "KV"). 

Hollerith input variable; ■ .''si '' . ■■r-!- 

= blank yields no effect 

= S when following item, D(5), is to be explicit t^^ 

= A when following item, D(5), is to be arctan t^ (in 

degrees) 

Hollerith input variable; 

= blank yields no effect 

= S when following item, D(6), is to be explicit t^ 

= A when following item, D(6), is to be arctan 

(in degrees) 

Hollerith input variable; segment shape (including NULL, 
in which case this segment is essentially deleted, and no 
further parameters are required) 

Hollerith input variable; segment type 

Value of 0 at the initial control point location for each 
arc (at the current x- station) 

Value of 0 at the final control point location for each arc 

Hollerith input; any comments 

x-component of surface unit normal 

y-component of surface unit normal 

z-component of surface unit normal 
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UTHETI(J) 


UTHET2(J) 

v(m) 

VTITLE(I) 

VX(M) 

VXX(M) 

W(I,J) 


WX(I,J) 


WXX(I,J) 

XCSMSl(KK) 

XCSMS2(KK) 

XDEL 

XEND 

XGO 

XNOW 

Y1(J) 

YIX(J) 

YIXX(J) 

Y2(J) 


Initial use 0 for each arc (as affected hy intersections 
and fillets) 

Final use 0 for each arc 

Current (latest x-station) computed value of each BL model 

I = 1 to 15; alphaniimeric vehicle or run title 

Current computed slope (dv/dx) of each BL model 

2 2 

Current computed derivative (d v/dx ) of each BL model, 

I = 1 to 4; defining mathematical parameters for each CS 
arc at a given station: 

R^(I = 1), 0^(1 = 2), A^(I =3), B^(I = 4) . 

I = 1 to 5; for I = 1 to 4, 

WX(I,J) = d(W(l,j))/dx 

WX( 55 J) = dr/dx for internal computations only 
I = 1 to 4; d(WX(l,j))/dx 

Starting x-station of the current cross section model 

Ending x-station of the current cross section model 

Input variable; increment size in x, to establish output 
stations between XGO and XEND 

input variable; final x-station to be interrogated 

Input variable; initial x-station to be interrogated 

Current x-station (used in various exercising routines; 
e.g., MODEl, M0DE2, etc). 

y of initial point for each CS arc 

,dYl(j)/dx 

d^Yl(<l)/dx^ 

y of final point for each CS arc 
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Y2X(J) 

Y2XX(J). 

Y3(J) 

Y3X(J) 

Y3XX(J) 

y4(j) 
y4x(j) . 
y4xx(j) 

ZCL(I) 

ZCLX(l) 

ZCLXX(l) 


ZMAPNM 

Z1(J) 

ZIX(J) 

ZIXX(j) 

Z2(J) 

Z2X(J) 

Z2XX(J) 

Z3(J) 

Z3X(J) 

Z3XX(J) 


dY2(j)/dx 

d^Y2(j)/dx^ 

y of slope control point for forward (initial) 
end of each CS arc 

dY3(j)/c^ 

d^Y3(J)/dx^ 

,y of slope control point for aft (final) end of each CS arc 

dY4(j)/dx 

d^Y4(j)/dx^ 

I = 1 to 3; current value (z) of bottom center line 
(l = 1), top center line (l = 2), and mapaxis (l = 3) 

I = 1 to 3; current slope (dz/dx) of bottom center line 
(l = 1), top centerline (l = 2), and mapaxis (l = 3) 

2 2 

I. = I to 3; current second derivative (d z/dx ) of bottom 
centerline (l = l), top centerline (l =2), and mapaxis 

(I = 3) 

Name of mapaxis 

z of initial point for each CS arc 

dZl'(j)/dx 

d^Zl(j)/dx^ 

z of final point for each CS arc 

dZ2(j)/dx 

d^Z2(j)/dx^ 

z of slope control point for forward (initial) end of each 
CS arc 

dZ3(j)/dx 

d^Z3(J)/dx^ 
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Z*+(J) z of slope control point for aft (final) end of each CS arc 

Z4X(J) dZ4(j)/dx 

24XX(J) d^Z4(j)/dx^ 





SYMBOL LIST FOR STEIN 


AAA, BBB, CCC, 

DDD,-EEE, FFF;"^ ^ __ . - — ■»- 

AAAZ, BBBZ, CCCZ, 

DDDZ, EEEZ, FFFZ, 

AAAZZ, BBBZZ, 

CCCZZ, DDDZZ, 

EEEZZ, FFFZZ The coefficients of the conformal mappings, and their first 
and second derivatives with respect to z 


ACH Free stream Mach number 

AEENF Dimensioneil free- stream pressure (Note: dimensions must 

be consistent with choice of length scale; this is for 
the computation of aero- coefficients only. ) 

Currently not used - leave blank . 

AR(I,J) I = 1 to KCOMP, J = 1 to KEIECE(i); integrated surface 

area for each component and piece 


AREF 

ARINF 

ATTACK 

B(M) 

BHH(M) 

BHZ(M) 

bn(m) 

bzz(m) 

B2 , B2Z 

C(M,L), CH(M, 
L), CZ(M,L) 


Reference area for aerodynamic coefficients 

Dimensional free-stream density (see note for APINF). 

Currently not used, leave blank . 

Angle of attack (input in degrees) 

Radial position of the body in the mapped plane 

Second derivative of body radius with respect to 0 

Cross derivative of body radius with respect to P and 0 

Radial position of body in mapped plane, at Z + DZ 

Second derivative of body radius with respect to P , 

y position (in the physical plane) of the wing tip and 
its derivative with respect to z (Fig. 5) 

Radial position of shock L in the mapped plane eind its 
derivatives with respect to 9 and P (mapped coordinates ) 
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CC(M,L),' 

ccy(m,l), 

ccz(m,l) 

cftitl(i) 


Radial position of L wing type shock surface (CC(m, 1), ■■ 
= B(M) and CC(M,L> = C(M,L-1) (L = 2 ... LC +l)) and its, 
derivatives with respect to Y and Z 

I = 1 to 5; alphanumeric request for aerodynamic coeffi- 
cients (i.e.j CL, CD, CM, CN, and CA). 


CMPTTL(l) 

CN(M,L), 

CHN(M,L), 

CZN(M,L) 

CONE 


DX(L) . 
DY(I) 

DZ 

DZFAC 

DZGEOM 

DZWRIT 

ERR(J) 

GAMFR 

GAMIN 

GAMLO (N,M) 

H(N,M) 

HCX(N,I),^ 

HCZ(N,I) 


I = 1 to kcOMP, CMPTTL (KCOMP + l) = TOIL (total); ' ' 

alphanumeric title for each component (above.) 

Radial position of shock L in the mapped plane at Z + DZ 
and its derivatives with respect to . 0 and ^ (mapped • 
coordinate at Z + DZ) r, , 

Cone half angle (input in degrees); only used for sharp cone 
calculations ; ;‘ 

Mesh spacing in the radial direction, in region L. , ' 

Mesh spacing in the circumferential direction, in regioh I. 

Step size 

Factor multipling DZ computed from CFL stability condition 
(usually DZFAC = .7) 

Interval for geometry test 

Interval for printed output 

error generated in an iteration 

Equivalent ratio -of specific heats (y) for frozen flow 

Free stream ratio of specific heats [y) ' 

2 ’ ' 

Local value of a /(p/p) 

(a== speed of sound, p == pressure, p = density) 

Mapped space polar angle 

th 

The X and Z derivatives of the I cross flow surface 
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HFN(I,J) 

HF0(I,J) 

HHL(M) 

HIN(I,J) 

HI0(I,J) 

ho(m) 

hs(n,i), 

HSR(N,I), 

HSZ(N,I) 

HSN(N,I), 

HSRN(N,I), 

HSZN(N,I) 

HST 

hi(m) 
HIN(M) . 

I 


lAERD 


Same as HFO(l,J.) but at ciorrent station (see Fig. 8) 

I = 1 to KC0MP, J = 1 to KPIECE(I); final value of 9' 
..def,ininp:,.each_com'Ponent_and^p iece„„a t ure vious. st ation ________ 

Value of 0 in the mapped plane of entropy layer surface 
point M (Fig. 3) 

Same as.HIO(ljj) but at ciirrent station (see Fig. 8) 

I = 1 to KC0ME, J = 1 to mECE(l); initial -value of 0' 
defining each component and piece at previous station 

Cylindrical 0' at. mesh points on the body at Z (see Fig. 8) 

Circumferential position of cross flow surface I and its 
derivatives with respect to r and P at Z 

Circvimferential position of . cross flow surface I in the 
mapped plane and its derivative with respect to r and P 
at Z + DZ 

Free stream total enthalpy 

Metric factor h^ (spreading of streamlines) at Z 

Metric factor h^ at Z + DZ • \ 

Counter for regions in the circumferential direction; 1=1 
in the region adjacent to the bottom symmetry plane and 
I = IC is the region adjacent to the top symmetry plane. 

. I is also a counter for cross flow type surfaces (l = 1, 
bottom symmetry plane; I = IC + 1, top symmetry plane). 

Indicator: 

lAERD = 0: Integrated forces and moments on the body are 

not read, and are set to 0. (This would be used to start 
an aero-coefficient run) 

lAERD =1: Integrated forces and moments on the body are 

read. (This would be’ used to continue an aero- coefficient 
run) 
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lAERO 


IBLOUT 


IBUG 


IC 

ICASE 


ICF(K) 


IDIMEN 

IENT(M) 


Indicator: 

lAERO =0: No aero-coefficients to be computed 

lAERO =1: At least one aero-coefficient to be computed 

Output (tape) unit for streamline/boundary layer code and 
sonic boom code - set equal to 0 if no boundary layer 
inputs are to be computed. 

Output indicator - IBUG = 0: no intermediate output, 

-IBUG= = 1: ' for intermediate output 

Number of regions in the circumferential direction 

Indicator - ICASE = 1: Initial flow field data are not 

read but computed in the code (i.e., first run for sharp 
nose vehicles ) 

ICASE = 2: starting plane data will be read (i.e., first 

run for blunt nose body or continuation run) 

K = 1 to 5; indicates request and name location for each 
aerodynamic coefficient (K = 1 for CL, K = 2 for CD, K = 3 
for CM, K = 4 for CN, and K = 5 for CA) 

ICF(K) = -1: coefficient not requested 

ICF(K) = N > 0: coefficient requested and 

CFTITL(N) = proper alphanumeric coefficient name (CL, CD> 
etc) . 

(If ICF(3) = 4, then CM is to be computed and CFTITL(4) = 
'CM') 

Maximum number of regions in the I direction 

Indicator for entropy layer IENT(M) = 0; surface not 
detected yet at M, IENT(m) = 1: surface detected at M, 

IENT(m) = 2: surface collapsed to body at M. 
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lENTE 


IFCP(I,J) 


IGAS 

IHS 


III 


INCP(I,J) 


IPUWCH 

IREADO 

IREADl 

IREAD2 

IREAD3 

IREAD4 


Indicator lENTE = 0: no entropy layer to be detected, 

lENTE = 1: entropy layer to be detected. lENTE is set 

equal to 2 when arv, entropy layer is started. 

I = 1 to KCOMP, J =1 to KPIECE(I); final control point 
(in 9) for each component and piece (determined from QUICK 
modeling) 

Indicator; IGAS = 0: ideal gas; IGAS = 1: equilibrium; 

is set equal to 2 at Z = ZREEZ (-freezing- station ) 

Indicator;. . . 

IHS =" 0: metric factor not computed 

IHS = 1: , initial plane data read and computed 

IHS = -1: initial plane data not read, but initialized 

by code to the body radius at each mesh point and computed 

Indicator: 

III = 0: No component pieces were found between this 

Z and Z + DZ 

III =1: At least one component piece was found between 

this Z and Z + DZ 

I = 1 to KCOMP, J = 1 to KEIECE(i); initial control point 
(in 0) for each component and piece (determined from QUICK 
modeling) 

Output \init for starting plane data for next run 


Set to 5 
Read unit 
Read unit 
Read unit 
Read unit 


in data statement in INIT-read unit for read 
for control data 1 
for control data 2 
for starting plane data 
for QUICK intermediate data 
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ISHEEG(J). 


Indicators: J = 1 denotes the bottom symmetry plane, 


ISHOK(M,L) 

ISHTIP 

IWRIT 
IZ(I,J) 

JA 
K 
KA 

KCOMP"' 

ICNTCAL 

KPIECE(l) 


■ J = 2 wing plane, and J = 3 the top symmetry, plane. 

ISHEEG( J ) =• 0 no sharp leading- edge at the J plane . 

■t" v> 

ISHBEG(J) - 1 there is a sharp leading edge at the 
plane but the shock has not been detected yet. ISHBEG( J ) 
is set equal to 2 vdien the shock has been detected. 
ISHBEG(j) is set equal to 3 when the shock is in. 

Wing type shocks surface indication for shock L at M 
ISH0K(M,L) := 0: arbitrary surface 

ISH0K(M,L) = 1: shock point (detached) 

ISH0K(M,L) = 2: sharp leading edge shock point 

Indicator: ISHTIP = 0 no sharp leading edges; 

ISHTIP 0 sharp leading edges exist on the geometry 

Output unit for printed flow field data 

I = 1 to KCOMP, J = 1 to mECE(l); Indicator: 

= 0: Component piece is not present between this Z and 

Z + DZ 

= 2: Component piece is present between this Z and Z ,+ DZ 

and, thus, must be integrated over 

Maximum number of steps between printed output 

Step counter, K = 0 at starting plane for each run 

Maximum number of steps before punching output and stopping 

run .-i. •'i-.-' ';nc!'vC''-' i.Mi ^ 

Number of individual components for -which aero- coefficients 
are to be computed 

Number of consecutive calls to AEROCF from ARCONT; signifi- 
cant for initialization procedures 

1=1, KCOMP; see NP 
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LC 

LDIMEN 

LOOP 


M 

MC(I), 

NSHKI(I), 

NSHK2(I), 

MREG(I) 

MCIR 

MCI 

MDIMEN 

MDZ 

MSHKl(L), 

MSHK2(L) 

MSHOK(N,l)," 


Counter for regions in the radial direction; L = 1 is the 
region closest to the body, L=LC is the region closest 
to the bow shock. L is also a counter for wing type 
shocks (L = — 1 inner most and L = LC — bow^shock). Finally 
L is used as a counter for radial dividing surfaces (i.e., 

L = 1 => body and L = LC + 1 => bow shock . ) 

Number of regions in the radial direction 

Maximum number of regions in the L direction - 

Indicator: 

LOOP = 0: level one of the MacCormack scheme 

LOOP = 1: level two of the MacCormack scheme 

LOOP = 100: print one more station and stop 

Counter in the circumferential direction; M = 1 is the 
bottom symmetry plane and M = MC(lC) + MREG(lC) is the top 
symmetry plane 

Correspond to NC(L), MSHKI(l), MSHK2(L). NEEG(L) but for 
cross flow type surfaces 


Minimum number of points in the "M" direction in any 
region I (usually MCIR =5) 

N\imber of points in the "M" direction in region 1 = 1. 
Maximimi number of points in the "M" direction 
The value of M at which the minimum step size was found’ ■ ■ 
Values of M at end shock points of shock L (Fig. 4-) 

Crossflow shock surface indicator 
MSH0K(N,I) = 0: arbitrary surface 
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w 

NC(L) 

WCl 

MDIMEK 

KDZ 

NLOOK 


NP 


NREG.(L) 

NRUW 

NSOUT 


P(N,M) 

PFT(I,J,K) 

PHL(M) 


MSH0K(N,I) = 1: cross flow shock point 

MSH0K(N,I) = 2; for points at a sharp leading edge shock 

Coimter in the radial direction (Fig. 2); N = 1 is the body 
and N = NC(LC) + NREG(LC) is the bow shock 

Number of points in region L (radial direction) 

Number of points in the radial direction in region L = 1 

Maximum number of points in N direction 

The value of N at which the minimum step size was found 

Indicator: 

= 0: wing type shock is first detected in any circum- 

ferential region I. 

= 1: wing type shock is first detected in region I =1. 

= 2: wing type shock is first detected outside of region 

1 = 1 . 

Number of pieces or segments into which a given aerodynamic 
component is to be divided (stored in KPrECE(l), I = 1 to 
KCOMP) 

NREG(L) = NC(L-l) + NREG(L-I) (NREG(1) = O) 

Run number, used to order runs 

Number of specific values of z at which there is to be 
printed output (NSOUT ^ lO) 

j6n(p/p^) (where p is the pressure) 

I = 1 to KCOMP, J = 1 to KFIECE(I), K = 1, 2, 3 ; X, y, and 
, z components, respectively, of the integrated pressure 
force for each component and piece 

£n(p/p^) on the entropy layer surface (Fig. 3 .) /' 
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phln(m) 

PIN 

PMT (I,J,K) 

PN(N,M) 

P0(M) 

po(n,m) 

R(N,M) 

rhl(m) 

rhln(m) 

RQRI 

S(N,M) 

SFR 

SHL(M) 

SHLN(M) 

sn(n,m) 

S0(N,M) 

T(N,M) 

TIN 

TRY(J) 

U(N,M), V(N, 
M), W(N,M) 

uhl(m), VH'^j 

(m), WHL(M, 


£n(p/p_j^) on the entropy layer siirface (Fig. 3) at Z + DZ 

p^/PgL (fi'ee stream pressure/sea level pressure) 

I = 1 to KCOMP, J = 1 to KPIECE(i); Cartesian components 
of the integrated moments for each component and piece 

£n(p/p^) (where p is the pressure) at Z + DZ 

In ARCONT and AEROCF only, £n(p/p^) on the body at Z 

i,n(p/p^) (where p is the pressure) at Z - DZ 

Mapped space radial coordinate 

Radial position of the entropy layer surface (Fig. 3) 

Radial position of entropy layer surface in the mapped 
plane at M and Z + DZ (Fig. 3) 

Ratio of the freezing plane gas constant to its free stream 
value 

Entropy 

Reference entropy at the freezing plane 
Entropy on the entropy layer surface (Fig. 3) 

Entropy on the entropy layer surface (Fig. 3) at Z + DZ 

Entropy at Z + DZ 

Entropy at Z - DZ 

Local value of (pressure/derisity) 

T^/Tgi (free stream teraperatur e/sea level temperature) 

J guess in an iteration 
Cartesian velocity components 

Cartesian velocity components on the entropy layer surface 
(Fig. 3) 
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uhln(m), 

VHLN(M)V 

whln(m) 

uw(nm), vn(n, 
m); wn(n,m) 

UN0R(I,J) 

VIN 

VMO(i)^ 

X(NW,L) 

XO(M) 

XTIP, XTIPZ 

Y(MM,L) 

YD, YDZ, YB, 
YBZ 

YQ(M). 

z ■■ ■ ' ' ■ ' 

ZCOMP 

ZEND 


Cartesian velocity" coniponents on the entropy layer surface 
(Fig. 3) at Z ■+ DZ 

Cartesian velocity components at Z + DZ 

The three (UNOR(l,J), UN0R(2,J) and UN0R(3,J) ) Cartesian 

. - < ' - ' "til 

components of the xinit normal to the body at the J sharp 

leading edge (Fig. 6) 

Free stream velocity" 

I = 2, 3; y and z positions of line about -which moments are 
computed 

Computational plane coordinate (x(l,L) = 0 and X(NC(L),L) 

= a) (Fig. 2) . -C ■ V , 

Cartesian x' at mesh points on the body at Z (see Fig. 7) 

X position (in -the physical plane) of the wing tip and its 
derivative with respect to z (Fig. 5) 

Computational plane coordinate (Y(1, I) = 0 and Y(MC(i),I) 

= 1) (Fig, 2) ■ ' = 

Position and z derivatives in the physical space, in the 
symmetry plane, of the top and bottom of the body. These 
roles depend on whether the configuration is hi^ -wing or 
low wing (Fig. 5) 

Cartesian y- at "mesh' points -.on the- body at Z (see Fig. ?) 
Axial station 

z station immediately prior to start of sharp leading edge 

Last axial station to be computed before punching output and 
stopping 




ZFINL(I,J) 


ZFREEZ 

ZGEOMl 

ZGE0M2 
ZINIT(I,J) 
ZMADD, MDEL 

ZMAPl, ZMAP2 


ZN 

ZNADD, NDEL 
ZO 

ZSHRP 

ZSOUT(l.) 

ZSTART 

ZTIPS 


I = 1 to KC0MP, J = 1 to KHECE(i); final station (z) for 
each component and piece (Note: ZINIT and ZFINL may over- 

lap or coincide for different pieces of the sanie component, 
thus allowing for disjoint cross sectional members) 

Value of z at which the thennodynamics is to be converted 
from equilibrium to frozen. 

First ^ial station at which a "geometry test" will be 
printed 

Last station of geometry test 

I = 1 to KCOMP, J = 1 to KPIECE(i); initial station (z) 

MDEL points will be added at Z = ZMADD. In the circumfer- 
ential direction 

The conformal mappings are not used for Z g ZMAPl and they 
are fully developed for Z ^ ZMAP2. (ZMAPl = starting plane 
station for the first supersonic flow run and ZMAP2 = ZMAPl + 
a number' of nose radii, usually) 

Updated axial station ZN = Z + DZ 

NDEL points will be. added at Z = ZNADD. ' In the "radial" . 
direction 

Z (in ARCONT and AEROCF) 

z station immediately following start of sharp leading 
edge (« ZCOMP) 

Specific .values at z- >at which there is to be printed output 
I = 1 -♦ NSOUT (if NSOUT ^ 0 no values of ZSOUT are read or 
stored) 

Starting value of z for run 

Value of z at which wing tip surface (Fig. 7) is inserted 
(usually ZTIPS < ZWING). This surface is used to control 
the grid. 
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ZWING 


ZWRITl 

ZWRIT2 

ZIMSH, 

ZINSH, 


Axial station. at which wing starts (used in mappings) 

(Fig. 7) 

Axial station at which output is begun (ZWRITl s ZSTART 
usually) 

Last axial station at which output is printed 
(ZWKLT2 < ZEND usually) 

Z2MSH Same as ZINSH, Z2NSH but for cross flow shocks (See Fig. 7 ) 

Z2NSH A wing type shock will be looked for between z = Z1NSH(J ) 

and z = Z2NSH(J). After detection, ZINSH(j) is set to 

6 6 

1 x 10 and Z2 NSH(j) is set to -1 x 10 so that shock J 

is not found again. (See Fig. 7 ) 



SYMBOL LIST FOR STRMBL 


ACHINF 

Free stream mach number, read from data tape 

ATTACK 

Angle of attack, in degrees, read from tape 

DZ 

Current step size, Az 

DZO 

previous step size 

FNU . 

, Nondimen sional kinematic viscosity 

, . - ■ 

FNUr=‘ = "'Joo'/'^ref ( ) = dimensional quantity and 


'^ref 

GAMMA 

Free stream ratio of specific heats, read from tape 

hcut(is,icut) 

6' -location of each streamline at each cut for body 
normals 

hp(n,m) 

Angle from x' -aocis (see Fig. 8) to mesh points (0' in 
Fig. 8) 

hpo(n,m) 

HP(N,M) at previous data plane 

HZNP 

d6'/dz for the current streamline and data plane 

HZOP(IS) 

d0'/dz for each streamline at the previous data plane 

hi(m) 

Metric coefficient h^ at mesh points on the body 

HIS 

Metric coefficient h^ for the current streamline and data 
plane 

IC 

Number of regions in the circumferential direction* 

ICO 

IC at previous data plane 

ICUT 

Indicator of current pseudo-stream-surface normal cut, 
from 1 to NCUT 


■31 


*As in STEIN 


ICUTMX 

Largest ICUT currently in storage 

IDUMl, IDUM2, 

Not used 

IDUM3 


IIC(IICUT) 

Indicates vdiich ICUT (= IIC(lICUT)) is currently stored in 
location referred to by IICUT 

IICUT 

Index (between 1 and NIICUT) to dynamic storage locations 
for pseudo- stream-siirf ace data 

INPT 

Index/counter for points taken" along' "body surface normals, 
from 1 to NNFT ' 

IR 

Read unit for card input 

IRT 

Not currently in use 

IS 

Streamline index/counter, from 1 to NS 

ITP 

l/O unit for data tape input 

IW 

Write unit for printed output 

JCUT 

Output and pseudo-stream-surface (p-s-s) parameter, normals 
to body are taken and p-s-s data is output every JCUT data 
planes 

JS 

Output parameter, streamline flow variables are output 
every JS data planes 

LC 

Niimber of regions in the radial directions* 

LCO ■ 

LC at previous data plane 

M 

Circumferential mesh point counter, from 1 to MC(lC) + 
MREG(IC) 

MC(I) 

I = 1 to IC; number of points in region I (circumferential 
direction)* 

MCO(l ) 

MC(l) at previous data plane 


*As in STEIN 
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MREG(I) 

MREGO(I) 

N 

NC(L) 

NCO(L) 

NCUT 

NFLG(lNra,IS, 

IICUT) 

NIICUT 


NNPT 

NREG(L) 

NREGO(L)’ 

NS 

NUM(IICUT) 


p(n,m) 

PI 


*As in STEIN 


MREG(l) = MC(l-l) + MREG(l-l), MREG(l) = 0 * 

MREG(I) at previous data plane 

Radial mesh point counter, from 1 to NC(LC) + NEEG(LC) 

L = 1 to LC; number of points in region L* 

NC(L) at previous data plane 

Number of pseudo-stream-surface normal cuts 

Flag set to indicate whether a point on the normal for a , 
given streamline has been computed (= 1) or not (= -l) 

Number of cuts permitted to be in storage simultaneously 
(must be sufficiently large, now equal to 5> to prevent 
body normal from the K + NIICUT data plane from extending 
past the K data plane or vice versa) 

Number of poihts taken along body surface normal to 
establish data in pseudo-stream-surface < 

NREG(L) = NC(L-l) + NREG(L-l), NREG(l) =0* 

NREG(L) at previous data plane 

Number of streamlines to be traced (up to 50} 

Number of points successfully computed for the IICUT cut 
(when NUM(IICUT) = NS*NNPT, all points on all normals 
taken at the ICUT corresponding to this IICUT have been 
computed, and thus may be output and the storage locations 
used for the next cut) 

£n(p/p^) at mesh points (where p is pressure) 
n 
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PN0RM(INPT,IS, 

IICUT), 

p/Poo (where p is pressure) at each point along the 
normal to each streamline for each cut currently 
being stored 

PO(N,M) 

P(N,M) at previous data plane 

PS 

An(p/p^) for the current streamline and data plane 
(where p is pressxire) 

RP(N,M) 

Radial distance from mapaxis (B^ line) to mesh points 
(r' in Fig. 7) 

RPO(N,M) 

RP(N,M) at previous data plane 

S(N,M) 

entropy at mesh points 

SLNG(IS) 

Integrated arc length along each streamline 

SN0RM(INPT,IS, 

IICUT) 

Entropy stored the same as PNORM(lNPT, IS, IICUT) 

SO(N,M) 

S(N,M) at previous data plane 

sr(is) 

r for" each streamline 

SS 

entropy for the current streamline and data plane 

STHE(IS) 

9 for each streamline 

TESTA 

Angle of attack, in degrees, read from cards 

TESTG 

Free stream ratio of specific heats , read from cards 

TESTM 

Free stream mach number read from cards 

TESTZ 

Initial value of z, read from cards 

THEOP(IS) 

9' for each streamline 

U(N,M) 

x-velocity component at mesh points 

UNORM(lNFT, IS, 
IICUT) 

X- component of velicty stored the same as PNORM(lNPT, 
IS, IICUT) 



UMX 


UNY 

UWZ’— - 
UO(N,M) 

US 

V(N,M) ■; 

VN0RM(INPI, 
IS, II CUT) 

VO(W,M) 

VS 

W(K,M) 

WN0RM(INPT, 

ISjIICUT) 

WO(N,M) 

WS 

YCL(J) 

YCLZ(J) 

YCLZZ(J) 

Z 

zcut(icut) 

ZO 

ZSTAR 


X- component of body surface unit normal 
y- component of body surface unit normal 

— z- component .of„body„surf ace^ unit^normal 

U(N,M) at previous data plane 

x-component of velocity, for the current streamline and 
data plane 

y-velocity component at mesh points 

y- component of velocity stored the same as PNOIM( INPT , IS , 
IICUT) 

V(W,M) at previous data plane - 

y- component of velocity for the current streamline and data 
plane . 

z-velocity component at mesh points 

z-component of velocity stored the same as PNORM(lNPT,ISy 
IICUT) 

W(N,M) at previous data plane 

z-component of velocity for the current streamline and 
data plane 

J = 1 to 3; y-position of body bottom center line (J = l), 
body top centerline (J = 2), mapaxis or line (J = 3) 

dYCL(j)/dz 

d^YCL(j)/dz^ - . o . 

Current z , ' 

z-locations at which cuts for body normals were made 
' z at previous data plane 
Initial value of z, read from tape 
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SYMBOL LIST FOR BOOM , 


ACHINF 

Free stream Mach number, read from data tape 

ATTACK 

Angle of attack, in degrees, read from tape ' 

DZ 

Current step size, Az 

GAMMA 

Free stream ratio of specific heats, read from tape 

HP(N,M) 

Angle from x' -axis (see Fig. 8) to mesh points (9' in 
Fig. 8) 

IC 

Number of regions in the circumferential direction* 

IR 

Read \init for card input 

IRT 

Not currently in use 

ITP- 

l/O unit for data tape input 

IW 

Write unit for printed output 

JA 

Output parameter, data are computed and output every JA 
data planes 

KZBDEX, KZTDEX 
KZCDEX 

See symbol list for QUICK (not used here) , 

LC 

Number of regions in the radial direction* 

M 

Circumferential mesh point counter, frcm 1 to MC(lC) + 

mreg(ic)* 

MC(I) 

I = 1 to IC; number of points in region I (circumferential 
direction)* 

MREG(I) 

MREG(I) = MC(l-l) + MREG(l-l), MREG(l) =0* 

N 

Radial mesh point counter, from 1 to NC(LC) + NREG(LC)* 

NC(L) 

L = 1 to LC; number of points in region L (radial 
direction')* 

*As in STEIN 
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NHPTS 

nreg(l) 

P(N,M) 

RCYL 

RP(N,M) 

S(N,M) 

TESTA 

TESTG 

TESTM 

TESTZ 

U(N,M) 

V(N,M) 

W(N,M) 

YCL(J) 

YCLZ(j) 

YCLZZ(J) 

Z 

ZSTAR 


Number of circumferential points on the data cylinder at 
which values of the flow variables are to be determined 

NREG(L) = NC(L-l) + NREG(L-I) , NREG(l) = Qx- 

J^n(p/p ) at mesh points (where p is pressure) 

GO 

Radius of data cylinder 

Radial distance from mapaxis (B^ line) to mesh points 
(r' in Fig. 8) 

entropy at mesh points 

Angle of attack, in degrees, read from cards 

Free stream ratio of specific heats, read from cards 

Free stream Mach number, read from cards 

Initial value of z, read from cards 

x-velocity component at mesh points 

y-velocity component at mesh points 

z-velocity component at mesh points 

J = 1 to 3; y-position of body bottom center line (J = l), 
body top centerline (J =2), mapaxis or line (j = 3) 

dYCL(j)/dz 

d^YCL(j)/dz^ 

Current z 

Initial value of z, read from tape 
. 1-1 ■■■ 


*As in STEIN 
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Figure h - CROSS SECTIONS IN THE PHYSICAL AND COMPUTATIONAL SPACES 
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Figure 7 - SHUTTLE ORBITER TOP VIEW 
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Figure 8 - AERODYNAMIC COEFFICIENT COMPONENT DEFINITION 
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INPUT DATA FORMAT 


INPUT DATA FORMAT FOR QUICK" " " """ " ^ 

QUICK input may be divided into three basic blocks: . data input for 
(l) cross section modeling, (2) body line modeling, and (3) exercising the 
model. The first block may also be subdivided into (la) - a cross section 
library definition, and (lb) - an application of this library to construct 
the total vehicle. For another presentation of QUICK input see Appendix A. 

(l) - Cross Section Modeling 


(a) - Library 



Card Type 

Format 

Variable Names 

1 

15A4 

vtitle(i) (I = 1, 15) 

2 

12 

NCSM 

3 

2I2,6X,10A4 

KDUM, KARC, CTITLE(i) 
(I = 1, 10) 


(Note: There will be exactly NCSM cards of type 3 appearing 

together with the appropriate cards of type 4.) 

4 A8,I2,A4,2X, ARCNAM, JSEQ, ASHAPE, 

a4,4x,ai,a8, atype, aspec(i), pntnam(i), 

■ , 1X,A1,4a8 ASPEC(2), PNTNAM(2), pntnam(3), 

arcnm(i), ARCNM(2) 

(Note: There will be exactly KARC cards of type 4 per model, and 
they will be grouped together for a given model after a card of 
type 3. ) 

(b) - Application (Note: These cards appear after NCSM blocks of one 

card 3 and KARC ceird 4 ' s . ) 


5 


I2,8X,A8 


KNTCSM, zmapnm 



Card Type 
6 


Format 


Variable Names 


212, 6X, KDUM, MODEL, XCSMBI(KDUM), 

2F10.5 XCSMS2(KDUM) 

(Note: There will be exactly KNTCSM cards of type 6.) 

(2) -Body Line Modeling 

Card Type Format Variable Names 

1 ai,a8,hc, byorz, bname, 

ai,a8 AYORZ, ANAME 

(Note: There will be as many cards of type 1, followed by its cards 

of type 2 and 3j as there are body line models, and as many cards 
of type 1, alone, as there are aliased control point coordinates, 
plus one blank card to terminate modeling input . ) 

2 I2,1X,a4, KSEG, SSHAPE, 

3X,AU,2X, STYFE,. SDEF, 

A2, II IFREE 

(Note: There will be as many cards of type 2 and 3 as there are 

segments in a given body line, plus one card type 2 with KSEG' = -1, 
These cards are deleted when aliasing.) 

3 3F10.5, D(l), D(2), D(3), 

2(F9.4', D(4), SLPl, D(5), 

Al), F10.5 SLP2, D(6) 

(Note: If SSHAPE is NULL, this card type 3 is deleted; also see 

Note for card type 2.) 

(3) - Exercising the Model 


Card Type 

Format 

Variable Names 

1 

12, IX, 

MODE, NDERV, XGO, 


12, 5X, 

XEND, XDEL, HGO, 


6F10.5 

HEND, HDEL 


(Note: MODE = 0, or blank, terminates all input.) 



An example of the input deck for a simple sharp-nose cone (10° half-angle)' 
with afterbody follows in Fig. 9* Figure Ilf also shows the intermediate 
data deck for this geometry. 

INPUT DATA FORMAT FOR STEIN 

There are five separate data sets read by the STEIN code. They are 
read on different read units because they may be generated in different 
places (i.e.j some may be user-generated and others are generated by other 
codes). These data sets are shown in Fig. 10. 

Control data (O) is read for every run of STEIN. This data set is 


generated by the user 
in INIT ) . The data in 

and read in on unit IREADO (set in a data statement 
control data (O) are 

Card No. 

Format 

.Variable Names 

1 

1615 

IREADl, IREAD2, IREAD3, 
IREAD4, IWRIT, IPUNCH, 
ICASE, IBUG, MCIR, NRUN, 
KA, JA, NLOOK, NSOUT, 

IBLOUT, lAERO 

» 


Control Data (l) is read for every run of STEIN. This data set is 
generated by the user and read in on unit IREADl. Its data are 


Card No. 

Format 

Variable Names 

2 

5F10.5 

ZEND, ZWRITl, ZWRIT2, 
DZWRIT, DZFAC 

3 

6F10.5 

ZGEOMl, ZGE0M2, DZGEOM, 
ZWING, ZTIPS, ZFREEZ 

4 

2(F10.5„I5) 

ZNADD, NDEL, ZMADD, MDEL 

5 & 5-a 

8F10.5 

ZSOUT(I) (I = 1, NSOUT) 
(if NSOUT < 0 these cards 
are not read) 
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Control Data (2) is read for every rim of STEIN. This data set is 
generated by the user for the first run of a configuration (geometry and 
free stream conditions). This data set is output (on IPUNCH) by STEIN 
for continuation rims of the same configuration but can be modified by the 
user. These data are read in on IREAD2 and consist of: 


Card No. 

Format 

Variable Names 

6 

5E15.5 

ZINSH(I) (1=1, 5) 

7 

5E15.5 

Z2NSH(I)..,(I ,= L, 5) . . 

8 

5E15.5 

ZIMSH(I) (I = 1, 5) 

9 

5E15.5 

Z2MSH(I) (I = 1, 5) 

10 

2E15.5 

ZMAPl, ZMAP2 

11 

715 

'■ lENTE, IGAS, ISHTIP, 
ISHBEG(I) (I = 1, 3) 
IHS 

The 

following data are 

read if and only if lAERO / 0 


Format 

Variable Names 


8X,I1,1X, 
5(A2, 3X) 

lAERD, CFTITL(I) (l = 1, 5) 


5E15.6 

VM0(2), VM0(3), APINF, ARINF, 
AREF 


12 

KCOMP. 


''■'I2,3X,AU 

NP, CMPTTL(l) (Note: NP is stored 

mECE(l)) 


12, IX, 12, 

INCP(I,J), IFCP(I,J), ZINIT(I,J) 


2F10.U 

ZFINL(I,J) (I = 1, KCOMP; J = 1, 
NP = KPIECE(l)) 
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The following data are read if and only if lAERD 7 ^ 0 (set and used by 
code for continuation runs ) . 

~~"~'~~~~~ Format" ^~~~~"”——~—— -- Variable s-Names — 

6 EI 3.7 PFT(I,J,K), PMT(I,J,K), 

AR(I,J) (I = 1, KCOMP; 

J = 1, NP = KPIECE(I); 

K = 1, 3) 

Starting plane control data are read for every run of STEIN. These 
data are generated by another code** or the user for the first run of a 
configuration. It is output from STEIN for continuation runs of the same 
configuration. These data are read on IREAD3 and consist of: 


Card No. 

Format 

Variable Names 

12 

4I5 

LC, IC, NCI, MCI 

13 & 1^+ 

5E15.5 

ZSTART, ACH, GAMIN, ATTACK, 



CONE, HN, TIN 

15 

3E15.5 

GAMER, RQRI, SFR 


(Only read if IGAS = 2 i.e., the 
flow has been frozen in a previous 
run of STEIN.) 

The starting plane flow field data are read by STEIN only if ICASE ^ Ij 
since if ICASE = 1 the starting plane flow field data are computed in STEIN 
(vehicle having a sharp circular nose of half angle CONE with axis the 
same as the Z axis). This data set is generated by another code** or the 

^ ; V J . V V' C T 

user for the first STEIN run and is output by STEIN for continuation runs. 
These data are received on unit IREAD3 and consist of: 


**These data are output by the BLUNT body code used to compute the flow 


over the nose of blunt vehicles. 



Format 

Variable Names 


NC(L), MSHKl(L), MSHK2(L), 
NEEG(L) (L = 1, LC) ' 

4I5. 

MC(I), NSHKI(I), NSHK2(I), 
MREG(I) (I = 1, IC) . ■ 

8011 

ISHOK (M, L) (L = 1, LC) 

(M = 1, MC(IC) +MREG(IC)) 


MSH0K(N, I) (I = 1, IC + 1) 
(N = 1, NC(LC) + NREG(IC)) 

4E13.5 

bn(m), CN(M,1), CHN(M,1) 
CZN(M,1) (M = 1 , MC(IC) 

+ MREG(IC)) ' 

3E13.5 

CN(M,1), CHN(M,L), CZN(M,L) 
(L = 2, LC) (M = 1, MC(IC) 

+ MREG(IC)) 

3E13.5 

HSN(N,I), HSRN(N,I), 
HSZN(N,I) (I = 2, IC). 

(N = 1, NC(LC) 4- NREG(LC)) 

5E13.5 

vn(n,m), un(n,m), wn(n,m), 
pn(n,m), sn(n,m), (N = 1, 
NC(LC) + NREC(LC)) and 
(M = 1, MC(IC) + MREC(IC)) 

The following data is read if 

and only if lENTE = 2 (i.e., entropy 

layer points have been detected): 


Format 

Variable Names 

Soil 

ient(m) 

(M =1, MC (IC) + MREG (IC)) 

6 EI 3.5 

RHLN(M), PHLN(M), UHLN(M), 

vhln(m), whln(m) shln(m) 

(M = 1, MC (IC) + MREG(IC)) 
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The following data are read if and only if IMS > 0 (i.e., metric 
coefficient is being computed, and is not to be initialized by the 
code ) . 

6EI3.5 HIN(M) (M = 1, MC(IC) + MREG(IC)) 

The QUICK intermediate data set is, read by STEIN for every run and 
is output by the QUICK code. These data are read on unit IREAD4. Since 
the user need not interact with these data, they will not be described 
in detail here. 

INPUT DATA FORMAT FOR STRMBL 

STRMBL input consists of user input control data, geometry data in 
the form of the QUICK intermediate data deck, and a flow field data tape 
generated by STEIN upon request. All control input is from unit IR, set 
in subroutine INOUT. 


Card No. 

Format 

Variable 

Names 

I 

4F10.5 

TESTM, TESTA, 

TESTG, TESTZ 

2 

EI3.6 

FNU 


3 

315 

NS, JS, JCUT 



Since the user need not alter the QUICK intermediate data deck, and 
the flow field data tape cannot be altered by the user, neither of these 
inputs need be described in detail. Geometry input is from unit IR; flow 
field data input is from unit ITP, also set in subroutine INOUT. 



INHJT DATA FORMAT FOR BOOM 

BOOM input consists of user input control data, geometry data in the 
form of the QUICK intermediate data deck, and. a flow field data tape 
generated by STEIN upon request. All control inputs are from unit IR, set 
in subroutine INOUT. 

Card No. ' Format Variable Name 

1 4F10.5 TESTM, TESTA, TESTG, TESTZ 

2 FIO.4, 215 RCYL, NHPTS, JA. . 

Since the user need not alter the QUICK intermediate data deck, and . 
the flow field data tape cannot be altered by the user, neither. of these 
inputs need be described in detail. Geometry input is from unit IR; flow 
field data input is from unit ITP, also set in subroutine INOUT. 
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Figure 9 - SAMPLE INPUT DATA FOR QUICK 
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Figure 10 - STEIN INPUT 
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OUTPUT FORMATS 


OUrMf FOR ^ 

QUICK generates several modes of printed output, output suitable for 
external plotting .codes, and an intermediate data deck (the inathematical 
model) to be used as input to other codes using SUB-QUICK. 

The math model is output on unit ITAPE (set in QUICK - the main 
routine) , from subroutine GEMOUT. ’ ITAPE may, of course, correspond to the 
punch unit in which case a card deck will be generated that may easily be 
used (with SUB-QUICK) with any other code. This data set need not be, 
altered (configuration changes should be made in the initial QUICK input 
data which should then be ' rerun ' through QJICK, -thus generating a new math 
model), and as such, will not he described in detail. This data deck is 
also included in the printed output and may be seen in Fig. Ilf. 

QUICK prints several cross section and body line checks with every 
'run. Fig. 11a .shows a correlation check between the- cross section input 
data and the math model. Labels and names make this and -all printed out- 
put self-explanatory. Note that. the indices in parentheses correspond to 
the indices in the tables. Any misspelled names will show up as additional 
items in the component and/or control' point tables and thus are easily 
detected on the first pass. A blank is always loaded into the first posi- 
tion of the control point index table. . ' 

Figure 11b shows a check list menu for body line models, output 
strictly for user convenience. In modeling a vehicle, the user may first 
define the logical cross section library and its application (see input 
data description) with subsequent blank cards to terminate input (thus, 
initially no body line models would be defined) and by filling in this 
table he could ensure that all control points were defined, either as a 
separate model or as an alias. 
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The output shown in Fig. 11c provides an important cross reference 
between the control point coordinates and the body line models (the 
indices in the parentheses) which define them. Model numbers are re- 
peated because aliasing was used. The left hand sequential index bears a 
direct relation to the control point index table in Fig. 11a. Each con- 
trol point has two coordinates viiich must be defined (y = f(x), z = g(x)), 
and in Fig. lie, the index for a particular control point's (n in Fig. 11a) 
z definition is m = 2n and for its y definition, m = 2n-l. Any control 
point coordinates that were not defined will have a zero (o) in the paren- 
theses, thus providing a quick check for complete definition. The first 
two blanks correspond to the initial blank in the control point index 
table of Fig. 11a. 

The output shown in Figs, lid and e provides a correlation check 
between the body line input data and the math model. The index in paren- 
theses represents the shape of that segment, a negative value indicating 
that a line between the initial and final points of that segment has a 
negative slope. The output of Fig. lie is completely annotated. In the 
column marked GAP, if two consecutive segments were not continuous in 
either x or v (v standing for y or z) the symbols X* or Y* would appear, 
accordingly. The last two lines in Fig. lie are generated in GEMOUT, and 
indicate that a successful check was performed to ensure that all control 
points are defined throughout the range of the cross section models in 
which they are to be used. 

Figure Ilf shows a listing of the math model. Figure 12 gives an 
example of the output, generated at user request only, from MODEl. The 
first line is an echo of the user's input which requested this exercising 
of the math model (MODE, NDERV, etc. ... see input data description). 

This line appears at the start of each piece of user requested output. 
INXBLM is the body .line model number, INXBLS is the segment number, and 
V represents y or z (VX = dV/dx, etc.). If MODE = - 1, no printed output 
will be generated, but the following output will be written on unit IPLOT: 



Line 


Variables 


Format 


1 


lAll, lANDV 

215 

2 


NXPTS, KNTBLM 

^ ^ 215 

3 


XNOW 

FIO.4 

4 


V(I), VX(I), 

3F10.4 



VXX(I) (I = 1, 




KNTBLM)- 



Blocks of lines 3 and 4 are repeated NXPTS times and line 4 is repeated . 
KNTBLM times for each line 3> • 

Figure 13 shows an example, of user requested output from M0DE2. The 
use of a "G" suffix (THETAG, RAD-G, ZGCORD) denotes variables referenced 
to the "geometric" coordinate system; i.e., the x .(not x') axis which does 
not include the shifting due to the mapaxis (most often the FRL), Note 
that in general, ZGCORD = Z-CORD + ZCL(3). Here, since z of the mapaxis . 
(ZCL(3)) is zero, ZGCORD = Z-CORD. Variables without the "G" are of 
course referenced to the mapaxis. H is used to represent 0', so RH, RX, 
RXH, and RXX are the first and second derivatives of the radius R with 
respect to 9' and x. • All labels with "SUB" indicate derivatives formed 
numerically in SLOPE. Where "SUB" appears together with "D" the variables 
shown are the differences between the analytically formed and numerically 
formed derivatives. Plotting output from M0DE2 (MODE = - 2) is in the 
following form; 


Line 

Variables 

Format 

. . J 

1 

lAMD, lANDV 

215 

2 

NXPTS, NHPTS 

215 

3 

XNOW 

FIO.4 

4. 

YP, ZPG, KNOW*, 
RPX*, RPH*, RPXX**, 
RPXH** 

7F10.4 


*written if and only if lANDV 2 1 

**written if and only if lANDV 2 2 
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Lines 1 and 2 are output once per call to M0DE2; line 3 is output NXPTS 
times per call; line 4 is output NHPTS times for each line 3* Line 4 out- 
put .consists of y, X, e, and r^g,. 

Output from M0DE3 is shown in Fig. l4. ZBCL, ZTCL, and ZMAP are 

ZCL(1)j ZCL(2), and ZCL(3)} respectively. J is an index reference for 

each arc, hut it may not appear sequentially since the arcs will he listed 

in increasing 0' after all intersections and fillets have been computed 

and inserted in their proper location. If J is positive the arc is in 

(IN(J) = l); if J is negative the arc is not in (lN(j) = - l) - this 

occurs, for exajiiple, when a growing piece is still completely contained by 

the basic skin or a fillet was unable to be inserted. U/THETAl and 

U/THETA2 are the theta limits of the arc, UTHETl( J ) and UTHET2(j) if 

J > 0, THETAl(j) and THETA2(J ) (original definition theta limits- unaffected 

by intersections or fillets) if J < 0. RO, HO, AA, and BB are curve para- 
2 2 

meters R^, ^ i B . The second portion of M0DE3 output is a cross- 

sectional interrogation in the neighborhood of each control point; labels 
are self-explanatory. 

Plotting output, for MODE = - 3 is generated in subroutine MODEl 
(multiple body line traces may be used to create plan and profile views). 
Output format is the seirae as for MODE = - 1 except for line 4 which will 
consist of just V(l), I = 1, KNTBLM (no VX(l) or VXX(l)). 

M0DE4 output is shown in Fig. 15. Labels are the same as those used 
in the output of M0DE2. 

Output from M0DE5 may be seen in Fig. I6. NORM-X, NORM-Y, and 
NORM-Z are the x, y, and z components of the unit normal to the body sur- 
face at the X, r' , 0' location indicated. 

There is also a mode of output for MODE = 6, but no separate sub- 
routine is involved. When MODE = 6 is specified, GEMCHK exercises modes 
1, 2, and 3 at x-stations near the limits of each cross secticn model. 

For plotting purposes, if MODE = - 6, GEMCHK exercises modes -2 and -7 
at these same stations. 
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M0DE7 output is for graphical purposes only. Output is again on unit 
IPLOT, and is in the form of cross-sectional cuts which show all arcs over 
their entire definition range (THETAl to. THETA2) rather than their limited- 
use range (UTHETl to UTHET2). For MODE - - 1 , output is in the following 
format : 


Line 

■Variables 

Format 

1 

lAMD, lANDV 

215 

2 ■ 

NXPTS, NHPTS 

215 

3 

KARC, KNTARC 

215 

4 

KNOW 

F10.5 

5 

Y, ZG, HNO^IR*, RX*, 

TT'10.5 


RH*, RXX**, RXH** 

Lines 1 and 2 are written once per call to MODE?. Lines 3 and 4 are 
written NXPTS times per call. Line 5 is written KARC*NHPTS times for 
each write of lines 3 and 4. NHPTS is the number of points on each arc, 
KNTARC is the total number of arcs at the current station, and KARO is the 
number of arcs minus any fillets that were unable to be defined at this 
station (and also the number of arcs output from this mode for plotting 
purposes ) . 

OUTHJT FORMAT FOR STEIN' 

STEIN generates three types of output. On unit IPUNCH STEIN will 
output (only if IPUNCH > 0 ) starting plane data to continue a run. This 
output is generated at Z = ZEND or at K = KA (i.e., the final axial station 
or step of a run). 


*written if and only if lANDV s 1 

**written if and only if lANDV ^ 2 
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The second type of output from STEIN is on unit IBLOUT (if and only 
if IBLOUT > O) and is used as input for both BOOM and STIMBL. IBLOUT 
should usually correspond to a tape unit, since a great deal of output 
is to be expected. This output consists of body and shock position, the 
flow field variables, and the various region sizing and control parameters 
(IC, LC, MREG(I), etc.) at each computational step. The formats are not 
important as long as they are consistent with the input formats of STRMBL 
and BOOM, and since all the formats are consistent they need not be 
discussed further. 

The last type of output from STEIN is usually printed o unit IWRIT. 
The input data is printed as shown in Fig. I7. The flow field data at 
the first axial station (Z = ZSTART) is always printed as in Fig. 16. 

Where X & Y are the Cartesian coordinates of the mesh point, P is the 
pressure (p/p^,) U, V & W are the three Cartesian velocity components, S 
is the entropy, M is the total Mach number and MA is the axial component 
of the Mach number. 'This flow field data will be printed in this format 
at every axial station between ZWRITl and ZWRIT2 at an interval of DZWRIT; 
the maximimi number of steps between outputs is JA. Figure I8 shows. a 
"Geometry Test" of the body in the mapped space. Here Y is the circum- 
ferential position in the computational space, B is the body radius in 
the mapped space, BH and BZ are the body derivatives with respect to the 
polar angle and axial position in the mapped space and finally XX and YY 
are the Cartesian coordinates ‘in the physical space. Figure I9 shows 
the output format for the variables on the entropy layer surface. 

Aerodynamic coefficients are also written on unit IWRIT following 
the flow field output at each z-station. An example of the aero-coeffi- 
cient output follows in Fig. 20a and b. The first piece of output, 20a, 
is computed using a reference area which is the integrated surface area 
of a given component up to the current station. The second piece of 
output, 20b is computed with a user input reference area. Labels make 
the output self-explanatory but it is important to note that the input 
reference area must be in the same units as the geometry is model. 


/ 
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OUTPUT FORMAT FOR STRMBL 

Output from STRMBL is of two main types. The first of these is 
associated with the tracing of _ streamlines on the body, and consists of 
the location (0', 0, r, x, and y) of each streamline and the value of the 
flow variables (u, v, w, p, and S) at these locations in various data 
planes. Also included are the index, the integrated arc length, and the 
value of the metric coefficient h^^ for each streamline at the current 
z-station, see Fig. 21a. 

The second type of output from STRMBL corresponds to the development 
; n 

of the pseudo-stream-surfaces. Locations and values of flow variables 

and their derivatives are output at MPT points along the body normals 

originating from each of the previously traced streamlines at selected 

data planes. For each data plane (which, along with the 9' location of 

each streamline and the geometry model, establishes the origin points for 

the body normals) there are two blocks of output associated with each 

streamline. The first block gives the location of and flow variable values 

at the points equally distributed along the body normal. The second 

block gives the length along the normal,, the derivatives of the flow 

quantities in the normal direction (DUDN = du/dn, etc.) and the corapon- 

/ ■* /N 

ent of velocity in the normal direction (VELDTN = Q • n or Q • Q) at the 
same points; see Fig. 21b. 

OUTPUT FORMAT FOR BOOM 

Output from BOOM, see Fig. 22, is a simple presentation of flow 
variables (p, S, u, v, w) on the surface of the data cylinder of user 
specified raidus with centerline at x = y = 0 (the z axis, not the z' 
axis). HC is the angle 9 to the points on the cylinder, measured from 
the windward symmetry plane. 
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Figure 13 - QUICK OUTPUT FOR MODE = 2 
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Figixre lit - QUICK OUTPUT FOR MODE = 3 
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Figure 15 - QUICK OUTPUT FOR MODE = 4 
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Figure 15 (continued) 
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Figure: l6 - QUICK OUTPUT FOR MODE = 5 
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Figure l8 - GEOMETRY TEST OUTPUT 
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Figure 19 - BOW SHOCK AND ENTROPY LAYER SURFACE OUTPUT 
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Figure 20a - AERODYNAMIC COEFFICIENTS OUTPUT I 
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Figure 20b - AERODYNAMIC COEFFICIENTS OUTPUT 2 
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Figure 21a - STREAMLINE OUTPUT" FROM. STRMBL 
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Figure 21b - PSEUDO STREAM SURFACE OUTPUT 
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STORAGE REQUIRElylENTS AND COMPUTER TIME 


STORAGE REQUIREMENTS AND COMPUTER TIME FOR QUICK' 

Using the IBM G-compiler, QUICK requires approximately 128 K^q bytes 
of core to compile (« words), and 176 K^q bytes to execute (« 54Kg 

■words).. CDC requirements may somewhat exceed the figures in parentheses 

since CDC machines do not use half-word instructions and IM machines do. 

. * 

These core requirements are true with the code dimensioned to allow 
a maxim-um of: 

10 arcs pre-cross section (maximum value of J*) 

10 segments per body line model (maximum value of N-**-) 

10 cross-sectional models (maximum value of K*) 

25 body line models (maxim'um value of M*) 

Of course, these may be adjusted if required. 

QUICK run time varies greatly with the user requested output options. 
On the IBM 370/l68, a sample run for a simple 10° cone with afterbody, 
exercising modes 1, 2, 3, ^ and 5 at four x-stations each, nineteen ( 19 ) 
circumferential points per station in mode 2, and seven circumferential 
points per station in modes 4 and 5 required approximately 30 cpu seconds 
(of which, less than a third would be attributable to the initial defining 
and checking tasks). On a more complex vehicle, exercising only mode 2, 
assembly of the model and output of data for thirteen cross-sectional 
stations, using theta increments of one degree (l8l points), required 
approximately 20 cpu seconds. 


*Each dimensioned variable in QUICK is defined in the Symbol list for 

QUICK in terms of these integers, unless otherwise specified. 



\ 

\ 

STORAGE REgJIREMENTS AND COMHJTER TIME FOR STEIN 

The ^storage used in STEIN is divided, of coiifse, between logic and' 
variables. Using fixed dimensions at a maximum grid of 40 x 50 (which 
... could be- required for 'very* complex vehicle's) the core needed to store the 
variables is iSOK^^ bytes (on the IBM 370/l68). The core required for 
logic without overlay is 400 K^q. So that 580K^^ bytes of computer core 
is needed to run STEIN in this configuration. When STEIN is overlayed, 
the core required for the logic becomes 160 K^q bytes. And if the dimen- 
sions of the variables were made to vary with/the problem the expression 
for core required for this part of the code would be (NDIMEN x MDIMEN) 

X 17 + MDIMEN X 70 + NDIMEN x 4o + 50Kj^q where NDIMEN is the number of 
points in the radial direction and MDIMEN is the number of points in the 
circumferential direction. For simple geometries with small shock layers 
these can be as small as 10 x 10. 

Presently the code is dimensioned to allow a maximum of: 

4o grid points in the radial direction (maximum value of N*) 

50 grid points in the circumferential direction (maximum value of 

M»«-) 

4 regions in the radial direction (maximum value of L*) 

4 regions in the circumferential direction (maximum value of I*) 

The computer time required by STEIN depends in general upon length 
of vehicle and free stream condition. One of the longest running calcula 
tions was that of a shuttle orbiter flying at M^ = 10 and an angle of 
attack of 30°. This calculation took about 2 hours on the GDC 66OO. 

Some of the reasons for this running time are: 

(1) At large angle of attack the shock layer on top of the body 
becomes large (requiring 25 mesh points in the radial direc- 
tion for accuracy). These mesh points are also across the 

*Each dimensioned variable in STEIN, STRMBL and BOOM is defined in the 
appropriate symbol list in terms of these integers, unless otherwise 
specified. 
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shock layer on the bottom of the body which makes, the physical 
distance between mesh points small and caused DZ '(- stable 
marching step)- to -become very small. With this small value r 
of DZ it takes 3000 steps to compute the entire vehicle. 

(2) On blunt nose vehicles the body entropy is very large causing 
small Mach numbers on the body. As the local aocial Mach number 
approaches one, DZ approaches zero. On the forebody of blvmt 
nose vehicles , this condition exists causing the calculation to 
slow down there . 

The computer time required to compute the flow field about an H.R.A. 
configuration at M^^^ = 6 and a = 0, was about 1 hour of CDC 66OO time. The. 
same number of mesh points at each axial station were computed in this 
case and the ^Shuttle arbiter case but the step size DZ -was doubled because 
of the small angle of attack and the low body entropy. Finally, the time 
required to compute the flow field about a simple slab delta wing 
and Q. = 30*^) from the nose to 15 nose radii down stream was about 15 min. 

The computer time/mesh points depend significantly upon two para- 
meters: 

(1) Vehicle geometry (Shuttle arbiter or simple slab delta wing) 

(2) Gas model used in thermodynajnics (ideal gas or chemical 
equilibrium) 

There is also a slight dependence on the number of imbedded shocks 
in the flow field, but this comparison is hard to make since one cannot 
run the same vehicle with and without imbedded shocks. 

STORAGE REQUIREMENTS AND COMRJTING TIMES FOR STRMBL 

With the IBM 370/l68 H-compiler, STRMBL requires roughly 2U0 K^q 
bytes of core to compile (;« words), and approximately 35 ^Kj^q bytes 

(w 131 Kq words) to execute. 
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Approximately eight cpu minutes were required to run STRMBL on the 
IBM 370/168 for an 89B shuttle calculation of about 225 computational 
steps (from Z = 50 to Z = 790; this piece of the flow field computation 
required approximately 22 cpu minutes.) 

STORAGE REgJIREMENTS AND COMPUTING TIMES FOR BOOM 

BOOM requires (for the IBM G-compiler) approximately 122K^^ bytes to, 
compile (w 37K0 words), and 190 Kj^q bytes to execute (« 60Kg words). 

In the same shuttle calculation as above, BOOM required about 3*6 cpu 


minutes . 



. PART 2 FROGRAJyiMER-ORIENTED DOCIJMENTATION ' 

OVERALL FLOW OF LOGIC 

QUICK consists of three basic sets of routines with distinct functions. 

The first of these reads the input data and begins to assemble 
the mathematical model - this is the defining portion of QUICK. 
The second set of routines perform some logical checking of 
the math model, and correlates it to the input data - this is 
the checking portion of QUICK. Included in this set is a rou- 
tine which reads user requests to exercise the math model, and 
calls upon the third and remaining portion of QUICK - the inter- 
rogating or exercising section, called SUB-QUICK in this report; 
see Fig. 23. 

STEIN utilizes a finite difference marching technique, so that given 
the flow field at one social station z the code computes the 
flow field at; z = z + Dz. This process is repeated until the 
desired station is reached. Figure 24 shows a flow chart of 
the overall logic used in STEIN. 

STRMBL performs two basic functions in two nearly independent steps. 

The first step reads all of the flow field data planes from 
tape and traces streamlines for the length of the vehicle in 
this run. Flow variables are evaluated and output along these 
streamlines. The link with the second step is the establishing 
of the cutting planes at which body surface normals will be 
taken to determine the pseudo-stream- surfaces (p-s-s). The 
data tape is rewound and control transferred to the second 
portion of the code which, reading through the entire data 
tape a second time, uses SUB-QUICK to establish the body normals 
and then evaluate the flow variables and their derivatives in 
the constructed p-s-s. An end-of-file (EOF) mark on the tape 
terminates the Job. 
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simply reads through the same flow field data tape used by 
STRMBLj and interpolates for flow variables on the data cylinder 
every JA data planes. (JA is a user input.) An end-of-file 
(EOF) mark on the tape terminates the Job. 















Figure 2k - STEIN OVERALL LOGIC 
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. OVERLAY DESCRIPTION 

The only code we found it necessary to overlay was STEIN. It '.was .■ 
found that the core req.uirernents could be reduced by 50 ^ using a simple 
overlay. 

The, routines in the root segment (No. l) , (always in core) are: 

STEIN (main routine), TIPSUR, • UPDATE, CSGEOM, . BLGEOMy CSC^C, IMAP,. 
MAP, BODY, NINTER, MINTER, PRAN, RANK, GAS, MOLEH, MOLES, EXPAN, 
OBSHK, SHTEST, SHTIP, VDOTV, MDOTV, THELIM, CSMINT, CSCALC , CURVES, 
'■ ■” CSMSET, CSMCOE, CSMFLT ■ 

Segment 2: INIT, GEOMIN 

' 3 BOUND ' ' ' ■ ■ 

4 SHARP 

. ' -5 ■ FREEZ . 

6 . NMESH ■ . ... 

■ 7 ' ■ ENTRLA 

' '■ 8 . SHMOVE ... 

9 MMESH ,■ ■■■.;■■■ ,>• 

10 ' . OUTPUT 

11 BLOUT 

12 POINTS 

"l3 COEF' ' ' ■ ' 

l4 NSHOCK ' ' ' 

r’ ■ 15 MSHOCK- . - , 

16 , MREGIO 

17 CFL 

18 SHREIN, SHPEDG 

'19' ARCONT, AEROCF,' KAREN 

20 NREGIO, INTSEC ' .... 

21 MSURFA, MTEST 

22 NSURFA, NTEST 
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SUBROUTINE DESCRIPTIONS 



SUBROUTINE DESCRIPTION FOR. QUICK 

BLGEOM assigns body line model values and derivatives to control point 
coordinates. 

BLMCHK correlates aind checks the; input data deck and the indices for 
.. -the generated body line math models . 

BLMDEF defines body line models from the input data. 

BLMSET controls the determination values and first and second derivatives 

for all body line models at a given x-station. 

CSCALC ' computes radial position and derivatives for specified cross 
section model, arc, and 0 '. 

CSGEOM is the main subroutine in the SUB-QUICK (look-up or exercising) 
portion of the QUICK system. It is called to establish r' = 
f(0^,x). It calls appropriate subroutines to evaluate body line' 
values and construct cross section geometry at a given x-station. 
It is used for all geometry model interrogation. 

CSMCHK correlates and, checks the input data deck and the indices for 
the cross sectional math model. 

CSMCOE composes the equations which are to define the cross section 
geometry at a given station. 

CSMDEF logically defines the cross section models frcan. the input data. 

CSMFLT creates control point definitions to permit the insertion of 

a smooth fillet between cross sectional arcs. 

CSMINT locates user specified intersections between cross sectional arcs 
and adjusts their use-theta limits. 
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CSMSET 


CURVES 

DLOKUP 


DSETUP 


GEMCHK 

GEMOUT 


GEOMIN 

KRVDEF 

MDOTV 

MODEl 

MODE2 

MODES 


sets up the control point coordinate arrays used to define 
the cross section geometry at a specified x-station. 

calculates values and first and second derivatives for individual 
curve fits. 

is a simple dictionary look-up routine. It assigns an index to 
match an input name to a codeword list, but is not capable of 
adding new items to that list. 

is an adapting dictionary look-up routine. New items are added 
to a codeword list, an index (counter) is returned for the 
codeword, and an indicator (INEW) is set equal to 1 \dien a 
hew item is encountered. 

exercises the mathematical model at user request via MODEl, 

M0DE2, etc. 

outputs the math model generated by the defining portions of 
QUICK (this is referred to as the QUICK intermediate data 
deck) . Also ensures that all body lines required by a cross- 
sectional model are defined. for the range of that model. 

reads in the math model generated by the defining portion of 
QUICK and output by GEMOUT (the QUICK intermediate data deck), 

calculates coefficients for the various curve fits associated 
with body line math models. 

performs matrix multiplication of a vector. 

is called by GEMCHK to trace body line model values. 

is called by GEMCHK to create cross sectional cuts. 

is called by GEMCHK to examine the cross sectional modeling 
in the region about control points. Mode -3 plotting is trans- 
ferred to MODEl (multiple body line traces to create plan and 
profile views). 
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modeU 

MODE 5 
MODE? 

QUICK 

SLOPE 

THELIM 

VDOTV 


is called by GEMCHK to exercise subroutine SLOPE and examine 
the numerically formed derivatives at various x-stations along 
traces at a- constant value o:D.‘e - .- - ---- - . 

is called by GEMCHK to examine the surface unit normals. 

is called by GEMCHK to examine all defined arcs at a given 
x-station. This routine is used for plotting' purposes only. 

is the main routine. It sets the read and write units and 
‘‘controls the flow of the defining, checking, and exercising 
• portions of the QUICK system. 

forms a numerical estimate of the first derivatives of a 
supplied set of points. It is used as an independent check on 
computed QUICK derivatives. 

creates and controls use-theta arrays to establish continuity 
in the cross sectional model. 

computes a vector dot product. 



SUBROUTINE DESCRIPTION FOR STEIN 

"AEROCF performs the integration ,of pressure forces and moments on 
the hody]^ for aerodynamic coefficient calculations. 

ARCONT controls the integration of pressure forces and- moments on the 
boc3y for aerodynamic coefficient calculations. 

BLGEOM (This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 

BLMSET (This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 

BLOUT outputs the entire flow field on tape at every computational 
step, to be used by STRMBL and BOOM. 

BODY computes the position (B(m)) of the body in the mapped space and 

its derivatives (BH(m) and BZ(m) ). The body is defined in the, 
physical space, in the routine BODY an iterative procedure is 
used to find the position of the body in the mapped space, and 
then BH(M) and BZ(M) are computed analytically. 

BOUND computes the position and derivatives of all boundaries of the 
computational space (CC(M,L), CCY(M,L), CCZ(M,L), HCZ(N,I) and 
HCX(N,I) ) from their positions in the mapped space. 

CFL computes the step size DZ that satisfies the Courant-Friedrichs- 

Lewy criterion for stability. It is called from the main rou- 
tine once per step. 

COEF computes the coefficients used in the conformal mappings and 

their derivatives. The positions of the top, bottom, and wing 
tip are transferred to COEF through common. These geometry 
variables are used to compute the coefficients of the mapping 
which are then stored in common. 

CSCALC (This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 
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CSGEOM 

CSMCOE 

CSMFLT 

CSMINT 

CSMSET 

CURVES 

ENTBLA 


EXPAN 


(This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines . ) 

(This routine is used both in- STEIN and QUICK, it is described 
in the section on QLJICK routines.) 

(This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 

(This routine is used both in STEIN and QUICK, it is described 
1 in the section on QUICK routines. ) 

(This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 

(This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines . ) 

is used to compute, detect, and collapse the entropy layer 
surface. It is called in each level of the MacCormack scheme 
(loop = 0 and LOOP =1). If lENTE is input as zero, control 
will return from ENTRLA immediately but if lENTE ^ 0 for the 
points on the entropy layer surface which have already been 
detected (IENT(M) = l) the position and dependent variables 
will be computed. When ENTRLA is called with LOOP = 1, after 
the dependent variables are computed, additional entropy layer 
points are looked for and all entropy layer points are tested 
to see which are to be collapsed (IENT(m) = 2) at the current 
station. 

computes the flow through a 2-D centered expansion corner. 
Given the upstream Mach number (XMl), GAML0(N,M) and the flow 
deflection (DELTA). EXPAN will compute the conditions after 
the expansion (pressure ratio P2QP1, temperature ratio T2QT1, 
Mach number XM2 and the slope (BETA) of the first expansion 
wave ) . 
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FREEZ 


GAS 


GEOMIN 


IMAP 


INIT 


INTSEC 


KAREN 


is called at a station Z = ZFREEZ when the thermodynamics of the 
flow field is in equilibrium. In FREEZ an equivalent "frozen 
state" is computed at each mesh point, IGAS is set to 2 so that 
the thermodynamics of the flow is frozen from that station oh. 
FREEZ is called, at most, once per vehicle. 

relates all the thermodynamic variables for ideal gas (iGAS = O), 
equilibrium air (iGAS = l) and frozen gas (iGAS =2). If 
IN = 1, P {JLn p/Po,) and S (entropy) are input; if IN = 2, P 
and H (enthalpy) are input; if IN = 3? S and H are input. GAS 
will compute GAMLO (N,M) and T(N,M) and then return if ICUT = 1. 
If lOUT ^ 1, GAS will compute the temperature (THE) and the 
variable P,S or H that is not input in addition to GAML0(N,M) 
and T(N,M) . 

(This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 

is the inverse mapping subroutine. It uses X and Y (physical 
Cartesian coordinates in the Z = constant plane) to compute R 
and THE (polar coordinates’ in the mapped space). The index I 
indicates which value of the coefficients (gotten in common) 

’'are to be used -- those at Z for 1 = 1, those at Z + DZ for 1 = 0. 

is used to initialize variables. In INIT all input data is read 
and then most variables are initialized. INIT is called only 
once per run. 

is called from NREGIO when two wing shock type shock points 
intersect. In INTSEC the conditions behind the resulting shock 
are computed. 

computes the area of the discrete triangular facets and sets up 
the unit normals used to integrate pressure forces on the body. 
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MAP 


MDOTV 

MINTER 

MMESH 

MOLEH 

MOLES 

MREGIO 

MSHOCK 

MSURFA 


is the mapping routine . It uses R and THE to compute X and Y 
(see description of IMAP) with the index I indicating at which 
value of Z the coefficients are to be used (as in IMAP). If 


ID = 0, X and Y are computed and control is returned. If 
ID = Ij the derivatives of the mapping, XR, YR, XZ, YZ, XH, 
YH‘(Xr, y^, y^, x^, y^) and. RX, RY, RZ, HX, HY, HZ, (r^, r 

r^, 6^5 9^, 0^) are also computed and returned in the argument 
list. In POINTS, for the body calculation; the second deriva- 
tives of the mapping are also needed, so that for ID = 2, 


y’ 


RXR, 

RYR, 

RZR, 

HXR, 

HYR, 

RXH, 

RYH, 

RZH, 

HXH, 

HYH, 

HZH, 

RXZ, 

RYZ, 

RZZ, 

HXZ, 

HZZ ( 

r , 

r , 

r , 

9 , 

6 > 

^ nl 

I" 






xr 

yr’ 

zr’ 

xr ’ 

yr’ 

X0’ 

y0’ 

ze 

®xe’ 


"■yP’ 


®x^’ 

Vp’ 


are 

computed and stored in 


common . 


(This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 


plays the same role as NINTER but for circumferential interpola- 
tion. I 

is called at Z = ZMADD to add MDEL points in the circumferential 
direction. These points will be divided proportionately between 
all the regions in the circumferential direction. 


uses curve fits of GAML0(N,M), T(N,M), S(N,M) and the tempera- 
ture as functions of P(£.n p/Pj^) and H (enthalpy) for air in 
equilibrium. 

uses an iteration to compute GAML0(N,M), T(N,M), H and tempera- 
ture (the) from P and S for air in equilibrium. 


shifts mesh points in the circumferential direction. There are 
no provisions for crossflow shocks intersecting. 

serves the same purpose as NSHOCK but for crossflow shocks. 

serves the same purpose as NSURFA but for crossflow shocks and 
surfaces . 
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MTEST 


serves the s^e purpose as NTEST but for crossflow shocks. . . 
Crossflow shock points started as infinitely weak shocks. 

NINTER- is. a general purpose interpolation routine. At some value of 

. M, NINTER interpolates- from an old mesh with NC(L) mesh points 

in LC regions onto a new mesh with NCN(L) • points in LCN regions. 
The positions of the old shocks are C(M,L) and those of the .new 
shocks are CN(M,L). . , 

NMESH. is called at Z- = ZNADD > to add MDEL points in the radial 

■ direction. These points will be divided proportionately betweai 
all regions in the radial direction. 

P^GIO' . >■ shifts mesh points in the:radial direction- as wing type shocks. 
approach each. other. ■■ When; two wing type shocks are .close 
enough to each other at some value of Y, they are intersected 
at that, point , the outer shock being considered the resulting 
shock .and the inner shock becoming an "arbitrary surface" at 
this point. When all the points on one shock intersect another, 
this shock is eliminated as a boundary. 

NSHOCK : computes the high pressure side of the wing type shocks, includ- 

: ... ling^the bow shock. NSHOCK is called from the main routine in 

each level of the MacCormack scheme.- After the interior points 
; - r;. have been ..computed in level one of the MacCormack scheme 
. -. NSHOCK uses the predicted values of the dependent variables on 
the low pressure side of the shock to integrate to a value of 
CZN(M,L). After level two of the MacCormack scheme the 
corrected values of the dependent variables on the low pressure 
side of the shock and CZN(M,L) compute in level one, are used 
to recompute the high pressure side of the shock. The bow 
shock is cpmputed only in level one since the flow on its low 
pressure side is constant. The position and derivatives (CH(M, 

L) and CZ(M,L)) of the wing shock type surfaces are also computed 
in NSHOCK. ’ 
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NSURFA 


NTEST 


ft:;'';.' 5 


OBSHK 


OUTPUT 


POINTS 


is used to rearrange the mesh when wing type shocks and wing 
shock surfaces are first inserted in the flow field. This 
routine is called after a shock point has been detected; in 
it the arbitrary surface is initialized. A new grid is defined 
and the dependent variables are interpolated. 

detects wing type shock points. If Z is not between ZlNSH(j) 
and Z2NSH( J ) , for some value of J, control is returned from NTEST. 
Once shock points are detected the initial jump conditions are 
gotten by extrapolating from either side and then CZN(M,L) and 
CHN(M,L) are computed. 

serves the same purpose as EXPAN but for a 2-D wedge compression 
. Both OBSHK and EXPAN are used in the sharp leading edge wing 
calculation. 

outputs on unit IWRIT the dependent and independent variables 
at each output station. The user specifies ZWRITl (initial 
output station), DZWRIT (output interval) and ZWRIT2 (last 
output station). The user can also specify NSOUT and ZSOUT 
for additional output. The maximum mmaber of steps between ■ 
output stations is JA and this routine will be called if execu- 
tion is terminated for any reason. When requested, aero- 
dynamic coefficients are also output from this routine. OUTPUT 
also writes (on unit IPUNCH) the starting plane data for the 
next run at Z = ZEND or K = KA (only if IPUNCH > O). 

computes all the dependent variables at interior points, body 
points, and on the low pressure side of all shock waves. For 
the portion of the internal boundaries that are not shocks 
the dependent variables are set equal across them in POINTS. 
POINTS is called from the main routine for each level of the 
MacCormack integration scheme. In POINTS the body second 
derivatives Btffl, BZZ, and BHZ are also computed. 
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PRAN •' 


RMK 


SHARP ; 


SHMOVE 


SHPEDG 


SHREEN 


SHTEST 


computes the flow through a Prandtl-Meyer centered expansion 
for equilibrium or ideal gas-. Given P(ln p/p^^) on either side 
of the expansion, the entropy (constant through the expansion) 
and the velocity in the plane of the fan (VNl) ■ PRAN computes 
the change in flow direction DXMJ. 

computes the flow through a shock. Given VWl (velocity normal 
to the shock), GAMl (the value of GAML0(N,M)), PI £n(p/p^), SI 
(entropy), Tl(p/p), and HI (enthalpy), all on the low pressure 
side of the shock, RANK computes these quantities on the high 
pressure side of the shock. 

computes the exact solution for the flow over a sharp circular 
cone at zero angle of attack (with half cone angle CONE) (for 
attached shocks). It also give an approximate solution for 
sharp cones at small angle of attack. SHARP is called once per 
run only if ICASE is input as 1. 

computes the positions and derivatives in the Z = constant plane 
of all shocks (crossflow and wing type). SHMOVE is called once 
per step from main. HN(N,M) is also computed here. 

computes the body unit normal components at a given fuselage 
station (x) on counterclockwise first (ILOHI = 1) or last 
(iLOHI =2) cross section arc ending or beginning with a 
control point at a specified angle (THE). 

iterates to find the exact location of the start of a sharp 
edge. Then it sets up a call to SHPEDG to establish the body 
normals . 

is used in the initial setup for starting a sharp leading edge 
wing. In SHTEST the mesh is adjusted to accommodate a sharp 
leading edge shock in the wing plane or top or bottom symmetry 
plane. 


SHTIP 


STEIN 

i i ■ 

THELIM 

TIPSUR 

UPDATE 

VDOTV 


calculates the flow variables behind a sharp leading edge wing. ■ 
In SHTIP, given the conditions in front of the sharp tip, the 
conditions behin^ the expansion^ or^compressi^n at the tip are_ 
computed. 

is the main program of this code. It is used for control mainly. 
In STEIN the geometry test is generated, some initialization is • 
performed, the marching loop is entered (i.e., ZN = Z + DZ) 
and finally, the routines that detect shocks or rearrange mesh 
points are called. 

(This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 

computes the position and derivatives (HSN(N,I), HSRN(N,I), and 
HSZN(N,I)) of the wing tip crossflow surface. 

is called once in each level of the McCormack scheme to "update" 
the dependent and independent variables. In UPDATE the symmetry 
conditions (U(N,l) = U(N,MC(lC) + MREG(lC)) = 0 and CH(1,L) = 
CH(MC(IC) + MREG(IC),L) = O) are also imposed. 

(This routine is used both in STEIN and QUICK, it is described 
in the section on QUICK routines.) 


95 



SUBROUTINE DESCRIPTION FOR STRMBL 


BLDEL 


BLGEOM 

BLMSET 

BRCKT 

BRCKTO 

CSCALC 

CSGEOM 

CSMCOE 

CSMFLT 

CSMINT 

CSMSET 


establishes the length of each line, in the direction of the 
body surface normal, which makes up the p-s-s. Currently this 
is an approximation for the boundary layer thickness on a flat 

plate 6 = 3 V-. .. = 5 * z/W Rez. 

(This routine is used both in STRMBL and QUICK, it is described 
in the section, on QUICK routines.) 

. • ■ ' « in. 

(This routine is used both in STRMBL and QUICK, it is described 
in the section on QUICK routines.) 

examines the distribution of mesh points in the current data 
plane to determine those points which will bracket a specified 
location. 


examines the distribution of mesh points in the previous data 
plane to determine those points which will bracket a specified 
location. 


(This routine is used both in STRMBL and 
in the section on QUICK routines.) • 

(This routine is used both in STRMBL and 
in the section on QUICK routines.) 

(This routine is used both in STRMBL and 
in the section on QUICK routines.) 

(This routine is used both in STRMBL and 
in the section on QUICK routines.) 

(This routine is used both in STRMBL and 
in the section on QUICK routines.) 

(This routine is used both in STRMBL and 
in the section on QUICK routines.) 


QUICK, it is described 
QUICK, it is described 
QUICK, it is described 
QUICK, it is described 


QUICK, it is described 
QUICK, it is described 
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CURVES 


DELTHE 

FLINE 

GEOMIN 

INCUT 

INTERH 

INTERR 

INTRHl 

INTR2D 

INTR3D 

LOCATE 

MAIN2 


(This routine is used both in STRMBL and QUICK, it is described 
in the section on QUICK routines.) 


controls the determination of flow variables on a given stream- 
line at the current station (Z), computes d6' /dz for the given 
streamline, integrates to find 0'g (circumferential location of 
the streamline) and (arc length measured along the stream- 
line), and determines r^ (radial position of the streamline). 


■ (e^, S^, and rg at Z + DZ).' 


is a simple function used for a line (where y = f(x)), determined 
from two distinct points, to calculate y* at a specific x*. 


(This routine is used both in STRMBL and QUICK, it is described 
in the section on QUICK routines . ) 

initializes all l/O units. 


performs a simple, second order interpolation in M (circumferen- 
tial direction) at a specified N. 

performs a simple, second order interpolation in N (radial 
direction) at a specified M. 

performs a simple, second order interpolation in M (circum- 
ferential direction) for variables only evaluated at the body . 
(a function of M only). 


performs a two dimensional, second order interpolation for 
quantities at a specified location. 


performs a three dimensional interpolation for any variable. 
The z-location of the point of interest must lie between the 
previous and current data planes. 


determines the location (z', r', 0') of a given point lying 
along the body surface normal taken at a specified z and 0'. 


is a subroutine, but acts as a second main program once STRMBL 
has established the z and '0' locations at which body surface 
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MDOTV 

WOUT 

sour 

STRMBL 

THELIM 

VDOTV 


normals are to be taken to establish the pseudo-stream-surfaces 
(p-s-s). The data tape is rewound just prior to entry into 
MAIN2, which then proceeds to search the flow field data, 
interpolating in three dimensions, and dynamically allocating . 
storage to find, store, and output all quantities of interest 
in the p-s-s. 

(This routine is used both in STRMBL and QUICK, it is described 
in the section on QUICK routines.) 

gives printed output of flow variables in the pseudo-stream- 
surfaces (p-s-s) and forms numerical derivatives of these 
variables in the p-s-s and outputs them. 

gives printed output of location and flow variable values for 
a given streamline. 

is the main routine. It reads data from cards and tape, calls 
the integrating and output routines, and sets up the stations 
at which the /Cuts will be taken for body surface normals to 
establish the pseudo- stream- svu:f aces. 

(This routine is used both in STRMBL and QUICK, it is described 
in the section on QUICK routines.) 

(This routine is used both in STRMBL and QUICK, it is described 
in the section on QUICK routines.) 
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SUBROUTINE DESCRIPTION FOR BOOM 

BLGEOM (This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

BLMSET (This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

BOOM is the main routine. It reads data from cards and tape, calls 

the appropriate interpolation routines, and outputs the data 
cylinder computed quantities. 

BRCKTl examines the distribution of mesh points to determine those 
points which will bracket a specified location. An INDEX is 
returned to indicate that the point was found in the field 
(index = O), inside the body (INDEX = l), or in the free stream 
outside the bow shock (INDEX =2). 

CSCALC (This routine is used both in BOOM arid QUICK, it is described 
in the section on QUICK routines.) 

CSGEOM (This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

CSMCOE (This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

CSMFLT (This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

CSMINT (This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

CSMSET (This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

CURVES (This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

FLINE is a simple function used for a line (where, y = f(x)), deter- 
mined from two distinct points, to calculate y* at a specific x*.. 
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. GEOMIN 


INOUT 


INTERH 


INTERR 


INTR2D 


MDOTV 


THEHTM 


VDOTV 


(This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

(This routine is used both in BOOM and STRMBL, it is described 
in the section on STRMBL routines.) 

performs a simple, second order interpolation in M (circumferen- 
tial direction) at a specified N. 

performs a simple, second order interpolation in N (radial 
direction) at a specified M. 

performs a two dimensional, second order interpolation for 
quantities at a specified location. 

(This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

(This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 

(This routine is used both in BOOM and QUICK, it is described 
in the section on QUICK routines.) 
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APPENDIX A 


A BRIEF CODE-ORIENTED USER'S GUIDE 
FOR THE QUICK GEOMETRY SYSTEM 



QUICK is a highly general geometry' package based on library controlled 
mathematical modeling of cross sectional arcs and body lines. The mathe- 
matical models for the cross sections and the defining lines. are taken 
together to provide a continuous analytic model of the surface geometryv 
Slopes, normals and all derivatives are therefore developed analytically. 

Of course, either discontinuous intersections or smooth fairings can be 
modeled and enforced in both the cross sections and the body lines. 

QUICK generally works in two basic coordinate systems (x, y, z) and 
(x, r, 0); see Figure Al. Data for modeling is input in Cartesian co- 
ordinates, while interrogations for exercising the models are performed 
in Cylindrical coordinates. Both of the coordinate systems are further 
subject to a translation in z. This is due to the necessary presence of 
a mapaocis, located in the symmetry plane, usually, corresponding to the 

position of maximijm half-breadth (y ); see Figure A2. The mapaxis is 
'■ mdjc 

necessary to fulfill one of the basic constraints of the QUICK approach, 
which is:' the radius (r) must be a single-valued function of the angle 
(6). Figure A2 (b) obviously does not meet this constraint, while 
\ Figure A2 (c), with a properly defined mapaxis, does. 

‘ During the discussion of the use of QUICK, several terms will appear 
frequently, and as such, will be defined here: 

' (l) . Cross section - standard definition; a planar cut through the 

vehicle normal to the FRL at a given x-station. 

(2) Cross sectional model - mathematical abstraction of a cross 
section, using simple curves to represent arcs between 
specified control points. 

(3) Control points - logically selected break or joining points 
between cross sectional arcs; initial and terminal points 
for idefining each arc. 



(4) Arc - a portion of one simple mathematical curve between two 
control points in cross section. 

( 5 ) Body lines - the defining lines of the vehicle geometry in plan’ 

and profile views; x-r\xnning control points given as - 

y^(x) and/or - z^(x). 

(6) Body line model - mathematical abstraction of a body line, using 
simple curves to represent segments between specified match 
points'. 

■ ■ ■' ' 

( 7 ) Match points - logically selected break or joining points 

between body line , segments ; initial and terminal points for 
defining each segment. 

(8) Se^ent - a portion of one simple mathematical curve between two 
match points of a body line model. 

( 9 ) Component - same as an arc; usually considered to be a nsimed 

portion of the vehicle geometry (e.g., a wing-upper-,^llipse 
may be component WNGUPELL). \ 

QUICK modeling is performed in terms of the basically independent 
logical cross section models and logical/mathematical body line mddels. 

The cross sections are defined purely in terms of the named componenji 

\ 

arcs and the named control points; see Figure A3 (a)-, which models the 
vehicle shown in Figure A2 (a). Body lines, corresponding to the named 
control points, are then defined mathematically for the length of the \ 
vehicle (or as long or short as is necessary); see Figure A3 (b). At a 
given x-station the body lines are interrogated to give values for the 
control points. These control point values are then used to create the 
required cross sectional arc models which are interrogated at a given 
value of 6. 

In cross section, a component arc is defined in terms of its control 
points, its shape, and its type. The arcs are considered to be ordered 
counter-clockwise (looking up the x-axis, i.e., in the negative x direction) 
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starting at the bottom of the vehicle (9 = - tt/2) and going to the top of 
the vehicle (6 = + n/2); see Figure A3 (a). A full complement of these 
arcs will define a cross sectional model, which is then given a specific 
range, in x, over which the model is applicable. The only exception, or 
extension, to the ordering rule is used to allow intersections between - • 

cross sectional arcs to be computed internal to the code. Components which 
may be considered to start in the body and grow out (such as a canopy or 
wing; see Figure A3 (a)) make use of AECNM, as defined later in Figure a 4, 
to. specify to the code the other arc sharing the intersection point. 

Such growing components are ordered as before except they appear after 
the last arc in the outer, basic skin. Fillets (see Table All and Figure 
Ah) are also ordered as before, but appear last as a group; i.e., all 
fillets follow both the basic skin and the growing adaptive pieces. 

The arc shapes available in cross section, along with their key 
words and equations follow in Table AI. 


TABLE AI - CROSS SECTION ARC SHAPES 


SHAPE 

KEYWORD 

EQUATION 

LINE 

LINE 

Ay + Bz + C = ' 0 

ELLIPSE 

ELU 

(y-y )^ . T 

(Concave to Origin) 


o + o = ± 

2 2 
A B 

ELLIPSEv 

(Convex to Origin) 

ELLO 

Same 


The line is defined exclusively in terms of its end points (control 
points); the ellipses also require a slope control point. 

The curve type controls the blending of the various arcs (or segments, 
since the cross sectional curves use the same group of curve types as' 
the body lines). In cross section, fore and aft are determined from the 
component ordering as mentioned before. A list of the curve types avail- 
able, their keywords, and their functions follow in Table All. 










TABLE All - CROSS SECTION AND BODY LINE CURVE TYPES 


(Blending Control) 


TYPE 

KEYWORD 

FUNCTION 

Piece 

PIECE 

Curve is defined as. a unit, with end 
points and slope control point if 
necessary. 

Aft-Link 

ALINK 

1 

*Curve being defined begins at the 
end of the previous curve and is 
tangent to it. 

Fore- Link 

FLINK 

1 

*Curve being defined ends at the 
beginning of the following curve and 
is tangent to it. 

Patch 

PATCH 

*Curve being defined begins at the 
end of the previous curve and ends 
at the beginning of the following 
curve and is .tangent to both of the 
adjoining curves. 

Fillet 

FILET 

End points and slopes of curve being 
defined are calculated from specified 
positions on the adjoining curves. 

**Null 

NULL 

Deletes an already existing segment. 

*In body line definition, "previous" and "following" are only relative, 
as the specific segments being linked or patched to are given as ., 
part of the data. 

**Available 

only in the modeling of body lines. 


Figure a 4, which follows, gives a card- by- card description of the data 
input format for cross sectional modeling. 

Consider, for an example, the simple forebody shown in Figure A5 (a). 
There are two basic cross sectional config\rrations corresponding to the 
initial purely conical section and the final section with flat sides. 
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One therefore selects the cross sections as shown in Figure A5 (h). The 
coding of the input data is shown in Figure a 6. Note that in the first 
model both ellipses are PIECE'S, while in the second model one ellipse 
is an FLINK and one is an ALINK. Also note the order in which the arcs 
are to be defined (JSEQ); for either of the ellipses to link to the line-,, 
the line must. first exist. Of course, depending upon the definition of 
the two slope control points, either or both of the ellipses could have 
been ' PIECE ' s . In the current setup, note that in model two the slope 
control points establish a slope for the center line points only, ^ the - 

slopes of the tangent points being established by the line. 

For a body line (a control point definition as a function of x), a 
segment is defined in terms of its match points, its shape, and its type, 
much . the seime as a cross sectional arc. The major difference between 
segment and arc definitions is that segment match points are numbers, 

i . - ■ 1 

establishing immediately the mathematical representation of the given • 
curve, while, as shown before, arc control points are, at the input stage, 
logical definitions only. Body lines may also be aliased to other body 
lines, when duplicate definitions are desired. The segments are consid- 
ered, to be ordered in the increasing x-direction over a range of applica- 
bility established by the match points. Segments are input in the order 
in which they are, to be defined and have an index to establish their 
x-direction ordering as opposed to the cross sectional arcs which are 
input in their order of appeareince (bottom to top) and have an index to 
establish their order of definition. This will be better understood after 
looking at Figure a 6 a little later and after having seen an example. 

A full complement of these segments (from one to the code's dimensional 
limits - these are presented later) will define a body line. 

■The segment shapes available are more numerous than are the arc 
shapes, and they follow, along with their key words and equations, in 
Table AIII. 


Il4 
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TABLE AIII - BODY LINE SEGMENT SHAPES 


_ BHAPE - . 

KEYWORD 

EQJATION . - 

Line 

LINE 

• Ax+By = O' 

X- Parabola 

XPAR 

Ax+By+y^ = 0 

y- Parabola 

YPAR 

. -n 2 ^ 

Ax+By+x =0 

Rotated X- Parabola 

RXPA 

Ax+By +Cxy+y^ = 0 

Rotated y-Parabola 

RYPA 

2 

Ax+By+Cxy+x = 0 

x-Ellipse 

ELLX 

2 2 

Ax+By +Cx +y .. = 0 

y-Ellipse 

ELLY 

2 2 

Ax+By+Cy +x =0 

Cubic 

CUBI(C) 

2 S 

Ax+By +Cx +x =0, 


The line is defined exclusively in terms of its endpoints; ,the x- and y- 
paraholas require, in addition, one slope to be specified and one to be 
left free; all other curves require two points, and two slopes (the slopes 
usually being established by means of a slope control point). T ■ 



The curve type controls the blending of the various segments ,. as 
for the cross sectional arcs. The list of curve types available for body, 
line segments, as well as arcs, along with their key words and functions; 
has already been tabulated in Table All. \ 

Following, in Figure a 6, is a card-by-card description of the data 
input format for body line modeling. A given segment is defined from an 
initial point as (x^, v^) to a final point (x^, v^) with an initial slope, 
tj^,- and a final slope, t^. Where applicable, t^^ and t^ are determined 
from a slope control point at (x^? ^ 2 )* letter "v" is used to 

represent y or z since either may currently be under definition. These 
cards follow the cross section data. 
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, Consider, for example, the same simple forebody that was used to 
demonstrate cross sectional modeling; Figure A5 (a). Looking back to our . 
cross sectional model, we see that we have defined a total of seven con- 
trol points (BDYBCL, BDYLTN, BDYLSCP, BDYUTN, EDYTCL, BDYUSCP and MAPAXIS). 
Each of these must now have y and z defined as a function of x. (The ma- 
paxis is constrained to the symmetry plane; i.e., y = 0.) Immediately 
following the cross section input data shown in Figure A7 one would input 
the body line data shown in Figure A8. Note that since tan (10°) = .176327 
the definitions for YBDYLTN and YBDYUTN are equivalent, and therefore 
could have been aliased.; Also note that in aliasing, only the model it- ’ ' 
self is important, and- thus one may alias ZBDYTCL with YBDYUTN.' Observe' 
that a negative reflection of a given body line requires a separate model. 

After readii.g the previous sections, a general approach to modeling 
any given configuration should begin to be apparent. One must first look 
at the general shapes involved in the cross sections, and determine how 
many unique cross section models are necessary to completely define the 
'vehicle. These cross sections must then be logically defined by choosing 
; the appropriate control points and . arcs as in Figure A3 (b) and Figure A5 ' 
(b), and by deciding upon each model's range of applicability, in x. The 
'coding of the input data for these cross sections can then be commenced. 
After this; one must carefully go through and define y(x) and z(x) for 
each control point. This completely defines the vehicle geometry. 

The code is currently dimensioned to allow 10 arcs per cross sectional 
model, 10 segments per body line model, 10 cross sectional models and 25 
■body line models. Of course, these may be adjusted if required. 

To exercise the geometry model, there are several modes of interro- 
gation built into QUICK. Following the blank card which terminates the 
body line modeling, one may insert a card of the format shown in Figure A9. 
A positive MODE produces printed output, a negative MODE produces a data 
file on unit IPLOT which may be used for plotting purposes. A blank card 
must follow these checkout requests to terminate the program. 

1 
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In the main routine, there ajpe five integer variables which control 
l/O operations. They are: 

_ ■_ ^ —.-IREAD-=— read-unit— 

iRlTE = write unit 

ICRITE = write linit for any error messages 
ITAPE = write/read unit for intermediate data file 
IPLOT = l/O unit for plotting data output from GEMCHK, 
MODEl,. M0DE2, etc. 


In addition, a reference piinch lonit (IPUNCH) is set equal to seven ( 7 ). in 
a data statement. This variable is used simply to' prevent improper l/O 
operations on the punch unit and is normally transparent to the user; how- ' 
ever, if the p\mch unit is not seven (7), then IPUNCH must be redefined' to 

the proper unit in QUICK and GEOMIN. ' “ 

■ ■ ■ ■ ■ , ■ 1 ■ ' 

1 

The intermediate data file is an interface between QUICK and SUB- 
QUICK. SUB-QUICK is a subset of QUICK' s subroutines which may be used iq • 
conjunction with any other code. In exercising QUICK, the intermediate 

. ’ i\ 

data file will be written on the unit corresponding to ITAPE. All neces- \ 
sary information is passed between the defining and checking subroutines 
and the interrogating subroutines of SUB-QUICK via common blocks when they ■ \ 

are used together; however, the intermediate data deck is both necessary \ 

arid sufficient for SUB-QUICK when exercising .it alone. A list of the 


routines in QUICK/SUB-QUICK follows: 

QUICK 

DSETUP 

DLOKUP 

CSMDEF 

CSMCHK 

BLMDEF 

/ " BLMCHK 

KRVDEF 

GEMOUT 

GEMCHK 

MODEI 

M0DE2 

M0DE3 
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mode4 

M0DE5 

MODE? 

SLOPE 

geomin' 

CSGEOM 

BLMSET 

CURVES 

CSMSET 

BLGEOM 

VDOTV 

MDOTV 

THELIM 

CSCALC 

CSMINT 

CSMFLT 


•SUB-QUICK 


To make use of SUB-QUICK, one must call two subroutines, the first 


being GEOMIN to read in the intermediate data deck, the second being 
CSGEOM for each point of interest. 


To read the data: 

CALL GEOMIN (ITAPE, IRITE, ICRITE, IREAD) 

Where: ITAPE = unit location of intermediate data deck 

for vehicle geometry 

IRITE = write unit 

ICRITE = write unit for any error messages 
IREAD = read unit (not currently used in SUB-QUICK) 
To interrogate at a point: 

CALL CSGEOM (X, H, R, RX, RH, RXX, RXH, NDERV) 

Where: X = x location 

H = theta location (-rr/2 ^ 0 ^ + n/2) 

R = radial distance from mapaxis to point on body surface 
corresponding to X and H. 

RX = dr/dx at this point 

RH = dr/d9 at this point 
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RXX = d^r/dx^ at this point 
2 

RXH = d r/dxd0 at this point 

NDERV = + N, __where_N_ is- the ordeh of derivative to.be calculated - 

+ N, previous call was to different location; must compute 

R and all temporary variables 

- N, previous call was to same point, thus derivatives 
may be computed without recomputing R or certain 
temporary variables 

The quantities X, H and NDERV are, of course, user specified, and the 

\ 

geometry code will return all other values. 

Two additional and more complex geometry modeling examples are 
included in Appendix A-A for the potential user's reference. 
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CNPYSCP 


.CNPYBTM 


(OTHER CONTROL POINTS 

AS be:fore) 



(OTHER CONTROL POINTS 
AS BEFORE) 



(a) 



Figure A3; - EXAMPLE OF QUICK DIVISION INTO 
CROSS SECTIONS AND BODY LINES 
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FORMAT 

SYMBOL 

DESCRIPTION 


Card 



1 



Col. 

1-80 

15A4 

VTITLE 

Vehicle or run title. 


Card 

_2 





Col. 

1-2 

12 

NCSM 

Number of distinct cross section models. 


Card 

(There will be 
of type U . ) 

exactly NCSM 

number of these cards appearing together with the appropriate cards 

Col. 

1-2 

12 

KDUM 

Running count of the current cross section model 

(from 1 to NCSM) 

Col. 

3-4 

12 

KARC 

Number of arcs in current cross section model. 


Col. 

11-50 

10A4 

CTITLE 

Title and/or descriptor of current cross section 

model . 

Card 

(There will be exactly KARC number of these cards per model; i.e. one for each arc, and they 
v/ill be grouped together for a given model after a card of type 3>) 

Col. 

1-8 

A8 

ARCNAM 

Arc or component name. 


Col. 

9-10 

12 

JSEQ 

Definition sequence (order in which the arcs are 

to be defined). 

Col. 

11-14 

a4 

ASHAPE 

Arc or component shape. 


Col. 

17-20 

a4 

ATYPE 

Arc or component type. 


Col. 

25 

A1 

ASPEC (1) 

= blank yields no effect . 



= Y when type is FILET, and only y - values are to be specified for 
that control point (z is computed on controlling component). 

= Z when type is FILET, and only z - values are to be specified for 
that control point (y is computed on controlling component). 

B to indicate that this control point is the bottom center line of 
the vehicle for this model (optional). 

- T to indicate that this control point is the top center line of the 
vehicle for this model (optional). 

Col. 26-33 A8 PNTNAM(I) Control point name for the beginning of this arc. 

Col. 35 A1 ASPEC(2) ■ Same as Col. 25, ASPEC (1). 

Figure At - DATA INPUT FORMAT FOR CROSS SECTION MODELING 


*l3T • 


Card 

_4 (Cont) 

FORMAT 

SYMBOL 

- Col. 

36-43 

A8 

PNTNAM(2) 

Col. 

46-53 

A8 

PNTNAM(3) 

Col. 

56-63 

a8 

ARCNM(I) 

Col. 

66-73 

A8 

;■ ARCNM(2) 




DESCRIPTION 


Control point name for the termination of this arc." - « .. 

Slope control point name for this arc when required, blank if not. 

if type is FILET; the name of the most aft component arc to which- 
the current arc's forward end is to be filleted. 

If type is other; the name of the most aft component ■ arc which, 
in case of intersection with: the current arc, 
is to update th'e forward end'bf the current arc 
and the aft end of the intersected arc. - . 

If type is, FILET;.-, the name of-; the most forward component arc to 

which' the current arc's aft end is to , be • filleted 

If ti^e is- other; the name of ’the most forward component _ arc which, 
in case of intersection with the current arc, is 
to update the aft end of the current arc and the 
; forward end of -the- intersected arc. 


Card type 3 
Col. 1-2 
Col. 11-18 


(Appears after NCSM blocks of one card type 3 and KARC cards of type 4.) - 

12 KNTCSM ■ Number of; cross section models to' define entire vehipler" 

A8 i' ZMAPNM' •' -' Name of mapaxis. - -i 

. •- . 't. 


Card type 6 
Col. 1-2'- 
Col. 3-4 ■ 


Col. 11-20 
Col. 21-30 


(There will be exactly KNTCSM number of these cards,) 


12 

KDUM 

12 

MODEL 


Flo. 5 XCSMSl 

F10.5 XCSMS2 


Running count of the current cross section model (from 1. to KNTCSM). 

Index corresponding to the already defined cross section models 
(between 1 and" NCSM ) . t ; J 

NOTE; KNTCSM may be larger than NCSM if a given model is, ‘used 
more than once. 

Starting x-station of the current cross section model. ■;/' 

Ending x-station of the current cross section model, 


Figure - DATA INPUT FORMAT FOR CROSS SECTION MODELING (Continued) 




100 ., 



. BDYTCL (BODY TOP CENTER LINE). 

-BDYUSCP (BODY UPPER SLOPE CONTROL POINT) i 


BOUPEL (BODY UPPER ELLIPSE) 



BDYLTN (BODY LOWER TANGENT) 


BDLOEL (BODY LOWER ELLIPSE) 


BDYLSCP (BODY LOWER SLOPE CONTROL POINT) 


BDYBCL (BODY BOTTOM CENTER LINE) 

BDYTCL 

BDYUSCP 

BDUPEL 


■BDYUTN (BODY UPPER TANGENT) 
■BDFLAT (BODY FLAT) 

•BDYLTN 

■BDLOEL 

BDYLSCP 

BDYBCL 


Figure A5 - SAMPLE FOREBODY WITH LOGICAL CROSS SECTION DEFINITION 
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Card 


FORilAT 

SYMBOL 

DESCRIPTION 

(Not 

e : There 

models 
blank 

will be 
, and as 
card to 

as many Cards of type 1, followed by its Cards of type 2 and 3> as there are body line 
many cards of type 1, alone, as there are aliased control point coordinates, plus one 
terminate modeling input.) 

Col. 

1 

A1 

BYORZ 

The letter Y or Z to indicate which data definition is to follow 
(a blank terminates all modeling input data). 

Col, 

2-9 

a8 

BWAl'-IE 

Body Line/Control Point name which is to be defined. 

Col. 

11 

A1 

AYORZ 

The letter Y or Z to indicate which definition is to be used when, 
aliasing (blank when not). 

Col, 

12-19 

a8 

ANAME 

Body Line/Control Point name to which BNAME is to be aliased, when 
applicable (blank when not). 

Col . 

31-70 

10A4 

TITLE 

Any comments , 


Card 2 (if not aliasing) 

(Note: There will be as many Cards of type 2 and 3 as there are segments in a given body line, plus one Card 

type 2 with KSEG = -1.) 


Col, 1-2 

12 

KSEG 

The order (in increasing x) in which this segment appears in this body 
line model. A KSEG = -1 (further arguments not required) terminates the 
data for a given body line (one Card 1). 

Col. U-7 

a4 

SSHAPE 

Segment shape (including NULL, in which case this segment is essentially 
deleted, and no further parameters are required). 

Col. 11-lU 

A4 

STYPE 

Segment type 

Col. 17-18 

A2 

SDEF 

Segment definition mode (currently, only two point, two slope/slope 
control point method is available - input "KV") 

Col. 19 

n 

IFREE 

Index of the datum quantity which is to be "free", i.e,, determined . 
by the code. IFREE ranges from 1 to 6 corresponding to Xj^,’ u^ , x , 
Vp, tj^, tg, as ordered. A line must have any one of these free; 
an X- or y(v)- parabola must have either 5 or 6' free; other curves 
should have IFREE =0. 


Figure 

a6 - DATA 

INPUT FORMAT FOR BODY LINE MODELING (Sheet 1 of 2) 
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F0K4AT 


SYMBOL 


DESGRIFTION 


Card 1 (if not aliasing) 

(see note for Card 2) 

(Note; If SSHA.PE is NULL, this 'card is deleted) 

Col. 1-10 F10.5 D (1) If type is PIECE, FLINK, this is x^. ■ 

If tjrpe is ALINK, PATCH, or FILET, this Is a floating point number i 
equal to KSEG of the segment from which and/or are to be determined. 

Col. 11-20 F10.5 ,D (2) If type is PIECE or FLINK, this is v^, ’ 

If type is ALINK, ’PATCH or FILET, this is a floating point number equal 
r . to KSEG of the segment from which t. is to be determined. 

Col, 21-30 F10.5 ;.D (3) If type is PIECE or ALINK, this is x^. 

If type is FLINK, PATCH, or FILET, this is a floating point number equal 

to KSEG of the segment from which x and/or v are to be determined!, 

2 d ‘ 'I ■ 

Col, 31-39 F9.U D (4) If type is PIECE or ALINK, this is Vg, jj 

If type is FLINK, PATCH, or FILET, this is a floating point number [equal 

* to KSEG of the segment from which tg is to be determined. ? 

Col, 40 A1 SLPl = blank yields no effect ! 

= S when following item, D (5), is to be explicit t^. 

= A when following item, I) (5), is to be arctan t^ (in degrees). 

Col, 41-49 F9,4 D (5) If SLPl is blank: 

If type is FILET, this is x^ (v^ and t^ are to be determined from the 
segment specified by D (l) and D (2)) 

. If type is other, this is x^ 

If SLPl is other than blank, see definition of SLPl, Col. 4o. 

Col. 50 A1 SLP2 = blank yields no effect 

= S when following item, D (6), is to be explicit t . 

= A- when following item, D (6), is to be arctan tg (in degrees). 

Col. 51-60 FIO .5 D (6) If SLP2 is blank: 

If type is FILET, this is X 2 (vg and t 2 to be determined from the 
segment specified by D (3) and D (4)). 

If type is other, this is 

If SLP2 is other than blank, see definition of SLP2. 

(i 

Figure a6 - DATA INPUT FORMAT FOR BODY LINE MODELING (Sheet 2 of 2) ; 
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FORMAT 


SYMBOL 


DESCRIPTION 


Col. 1-2 12 MODE = 0 (or blsmk), terminates all input. 

= ± 1, creates body line traces. 

= ± 2, creates cross sectional cuts. 

. ■ It 

= + 3> interrogates cross sections in neighborhood of control points!;. ' 

. I 

- 3j allows multiple body line traces to create plan and profile views. 

= + 4, comparison of analytic derivatives with numerically formed 
derivatives. 

= + 5j check of unit vectors normal to body s\rrface. 

= + 6, exercises modes 1, 2, and 3 at the limits of each cross 
sectional model. 

- 6, exercises modes - 2 and - 7 at the limits of each cross I 

sectional model. , ' 

i 

= - 7, (plotting mode only) creates cross sectional cuts, but includes all 
arcs in their entirety (including growing pieces still contained 
within the basic skin). 

Col. 4-5 12 NDERV = ± N, where N is the order of derivative to be calculated (N=0, 1 or 2). 

+ N, should always be used for these interrogations (means each call to 
a given location is new, thus the radius and all temporary variables 
must be computed). 

- N, should not be used for these interrogations (requires previous call 

to same location (x and 0), radius and certain temporary variables are 
not recomputed). | 


Col. 11-20 

• F10.5 

XGO 

Initial x-station to be interrogated. 


Col. 21-30 

FIO .5 

XEND 

Final x-station to be interrogated. 


Col. 31-40 

FIO .5 

XDEL 

Increment size in x, to establish outputs stations between XGO and XEND 

Col. 41-50 

Flo. 5 

HGO 

Initial value of theta (in degrees) to be interrogated; not 
for modes 1, 3. 

required 
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Final value of theta (in degrees) to be interrogated; not required : 
for modes 1, 3. j 

Col. 61-70 
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HDEL 

Increment size in degrees. to establish interrogation points 

between ‘ 


HGO and HEND; not required for modes 1, 3. 
Figure' A9 - DATA INPUT FORMAT FOR EXERCISING THE GEOMETRIC MODEL 
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APPENDIX B 


A BRIEF USER'S GUIDE 
TO THE 

THREE-DIMENSIONAL BLUNT BODY CODE' ( BLUNT.) 



BLUNT is a simple to use code which will accept the QUICK intermediate 
data deck to define a hlunt nose body and will supply a directly useable 
(data deck for the starting plane of STEIN. 

Here BLUNT 's input data will be described. There are three input 
data cards for BLUNT in addition to the QUICK INTERMEDIATE DATA DECK. 
This intermediate data deck is output from QUICK and the user need not 
get' involved in its details. 

Input : 

Card #1 NRUN, MONTH, MDAY, MYEAR, NA, MA, LA, KA, JA, LB, LE, IN, 
IGAS, IRESTRT 

Card #2 ACH, GAMMA, STAB, THEMAX, ELL, XO, ANGLE, ALPHA 
Card #3 PIN, TIN 
Card #4 NCSU, MCSU, IPUNCH 
QUICK INTERMEDIATE DATA DECK 
Formats ; 

All quantities on Card #1 are read in 15 format 

All quantities on Cards # 2 . and #3 are read in E10.4 format 

All quantities on Card #4 are read in 15 format 

Nomenclature : 

NRUN Run number 

MONTH Month 

MDAY Day , 

MYEAR Year 

NA Number of intervals in the r direction (maximum of 10) 

(Fig. Bl) 

\ 

MA Number of intervals in the 9 direction (maximum of lO) , 

. (Fig. Bl) 
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LA 


Number of intervals in the 0 direction (maximum' of 8 ) 

(Fig. Bl) 

KA The number of steps to be computed, after which the code will 

output initial data. Typically KA = 700 to reach steady state. 

JA The n\imber of steps between outputs before the steady state. 

LB Indicator for output quantities indicating convergence at every, 

step LB = 0 for no output. 

LE Geometry indicator: 

LE = 0 General geometry input (from "QUICK") 

LE ^ 0 Circular cross sections, geometry is nondimensionalized 
with respect to 'Che radius of curvature of the nose. 

LE = 1 Paraboloid cap 

LE = 2 Ellipsoid cap with a given axis ratio (ELL) and 
followed by a cone with half angle (ANGLE) 

IN Index not used 

IGAS Gas Indicator IGAS = 0 for perfect gas IGAS = 1 for air in 

equilibrium 

IRESTRT Restart indicator: 

I 

= 0 Blunt body is started with code supplied guess and outputs 
data on unit 8 for restarting blunt body code. 

= 1 BLUNT reads starting data from unit 8 and continues. 

‘l . 

ACH Free stream Mach number 

GAMMA Ratio of specific heats (Cp/Cv) in free stream 

STAB Stability factor for C.F.L. condition (DT = DTmin(STAB) ) . 

Typically STAB =1.2. 

Limit on 9 . Now computed in code but still in read statement. 


> 

1U7 


THEMAX 



ELL 

XO 

^GLE 

ALPHA 

S i 

;piN ^ 

;tin ■' . 
ijjcsu 

MCSU ^ 
IPUNCH . 


Used only when LE = 2. Axis ratio of ellipsoid, ELL = 1. for 
spherically caped cone. 

Location of center of coordinate system (Fig. Bl) (XO should be 
large enough so that initial data plane for supersonic flow 
calculation in supersonic) 

.Cone half angle for LE = 2 V ... 

Angle of attack 

Free stream pressure (p /Pot ) ’Jse only when IGAS = 1. 

oJj 

jr. * . 

Free streami temperature (T^/Tg-^^) used only when IGAS =1. 

Number of mesh points in the initial data plane in the r 
direction (Fig. B2). NCSU can be different from NA + 1. 

Number of niesh points in the 9 direction (Fig. B2) 

Output unit for initial data plane results. 




Figure B1 - COORDINATE SYSTEM 



.Figure B2 - ILil^riAL DATA PLAICE •’ 
1^9 ■ 
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